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I .  INTRODUCTION 


It  has  been  known  for  many  years  that  the  geomagnetic  field,  particularly 
under  nighttime  propagation  conditions,  causes  the  normal  modes  in  the  very 
low  frequency  (VLF)  and  low  frequency  (LF)  bands  to  be  hybrid  (ref.  1).  That 
is,  the  modes  are  neither  pure  transverse  magnetic  (TM)  nor  pure  transverse 
electric  (TE)  but  rather  a  mixture  of  the  two  polarizations.  Thus,  a  vertical 
electric  dipole  which  radiates  pure  TM  wave  will,  by  virture  of  the  radio  wave 
interaction  with  the  ionosphere,  be  a  source  of  TE  wave  as  well.  At  the 
ground  the  TE  field  is  normally  constrained,  because  of  relatively  high  ground 
conductivity,  to  be  very  much  less  than  the  vertical  electric  field  component. 
Nevertheless,  the  longitudinal  rf  magnetic  field  component  at  the  ground  can 
be  significant  relative  to  the  TM  component.  This  is  possible  since  the 
longitudinal  rf  magnetic  field  component  depends  upon  a  height  derivative  of 
the  TE  field  component.  The  longitudinal  magnetic  field  component  can,  in  the 
case  of  direction  finding  with  loops,  cause  appreciable  direction  of  arrival 
(DOA)  error;  especially  in  proximity  to  ranges  where  the  transverse  magnetic 
field  undergoes  a  modal  interference  null.  Though  capability  for  estimating 
DOA  errors  due  to  the  geomagnetic  field  has  existed  for  sometime,  the  need  for 
such  calculations  was  minimal.  New  VLF  and  LF  digital  receivers  under  Navy 
development  will  include  antenna  pattern  steering  algorithms  whose  performance 
depends  upon  DOA  properties  of  the  received  signal  (ref.  2).  This  makes 
timely  the  modification  of  VLF/LF  propagation  codes  for  the  purpose  of 
predicting  DOA  errors  associated  with  geomagnetic  field  effects. 

Included  in  the  present  report  is  a  modification  of  an  early  mode  conversion 
program  (ref.  3)  which  allows  for  error  calculations  of  the  DOA  as  a  function 
of  range  for  both  uniform  and  laterally  inhomogeneous  guides.  As  with  the 
earlier  program,  the  present  program  requires  waveguide  modal  constant  inputs 
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as  well  as  quantities  relating  to  excitation  factors  and  modal  polarizat i on . 
These  inputs  are  standard,  NPUNCH  ---  1,  output  of  the  waveguide  program 
"MODESRCH”  (ref.  4).  Principal  outputs  of  the  present  program  are  mode  sum 
plots  as  a  function  of  range  for  scaled  transverse  and  longitudinal  magnetic 
field  components  along,  with  a  DOA  error  curve.  Provision  is  also  made  for 
plotting  these  quantities  for  fixed  t  ransmi  ttei'- receiver  distance  as  a  func¬ 
tion  of  position  of  a  horizontal  inhomogene  i  t  v  such  as  the  day-night 
terminator.  This  option  is  useful  only  if  the  ground  conductivity  and  geomag¬ 
netic  i  i  e  1  d  magnitude  and  orientation  may  he  regarded  as  constant  over  the 
path.  It  is  emphasized  that  the  present  program  is  restricted  to  these  out¬ 
puts.  Should  the  usual  electric  field  components  be  desired,  the  user  is 
directed  to  an  update  of  the  program  of  reference  3  (which  can  be  obtained 
from  L.R.  Hitnev)  or  to  the  program  of  reference  5.  The  latter  program  can 
also  be  easily  modified  to  accommodate  DOA  capability. 

Reference  2  provides  results  of  ground  based  measurements  at  San  Diego  of  the 
variations  in  the  DOA  of  VLF/LF  signals  from  a  number  of  fixed  vertical  trans¬ 
mitters.  A  major  purpose  of  the  present  study  has  been  to  compare  computer 
result's  with  these  measurements.  Though  reference  2  provides  results  for  24- 
hour  periods,  only  the  all-night  and  all -day  cases  are  considered  here. 

In  sect  ion  11  the  origin  of  the  rt  longitudinal  magnetic  field  component  is 
briefly  discussed  and  relevant  formulas  are  summarized.  A  brief  description 
of  the  program  is  given  in  section  Til.  For  a  more  complete  discussion  of  the 
program  the  reader  is  referred  to  references  3  and  6.  Results  are  presented 
in  section  TV  and  conclusions  are  given  in  section  V.  The  Appendix  contains  a 
p rograii;  1  i  r;  t  i  ng . 
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II.  SUMMARY  OF  EQUATIONS 


In  the  following,  an  x,  y,  z  rectangular  coordinate  system  is  assumed  with  x-z 
the  plane  of  propagation  and  positive  z  directed  into  the  ionosphere. 
Invariance  in  y  is  assumed.  Waveguide  theory  applied  to  VLF/LF  yields  in  the 
asymptotic  limit  (refs.  7  and  8) 


_ K _ 

[ sin(x/a) ] 1 '/2 


exp( -iksin(#m)x) 


(1) 


where  K  is  a  constant,  m  is  a  mode  index  and, 
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=  y  component  at  the  rf  electric  field. 


field  intensity  (H^  =  Z where 
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=  earth's  radius . 

=  receiver  height. 
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=  modal  height  gain  for  y  component  of  the  rf  electric  field. 
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=  free  space  wave  number. 
=  modal  eigenangle. 
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For  laterally  Homogeneous  guides  A^  is  the  excitation  factor  which  depends 

upon  the  dipole  transmitter  altitude  and  orientation  as  well  as  upon  functions 
of  the  eigenangle.  For  laterally  inhomogeneous  guides  A^  depends  as  well  upon 

mode  conversion  coefficients.  Details  of  the  functional  dependencies  can  be 
found,  for  example,  in  references  3  and  6.  However,  those  details  are  of  no 
particular  significance  in  the  present  discussion.  What  is  of  consequence  is 
the  origin  of  e^  for  a  vertical  transmitter.  To  understand  this,  use  is  made 

of  the  consistency  relation  which  requires  that  a  downgoing  (d)  plane  wave  at 
the  ground  must  after  ground  and  ionospheric  reflections  return  to  its  same 
value  upon  returning  to  the  ground.  Mathematically  this  may  be  written  as: 


a-uejuej) 


hd  (o) 
ymv  ' 

d  ,  N 
,  e  (o) 
^  ymv  ' 


-  0 


where 


(2) 


R(fl)  - 


I,R|,  <*)  xRll  <*>' 
llR1(5)  ^(0) 


(3) 


is  a  plane  wave  reflection  coefficient  matrix  referenced  to  the  ground  and 
represents  the  reflection  from  everything  above  ground  level  with  vacuum  below 
it  and 

nR.i  (9)  0 

o  iRi(*>  (4) 


R(tf)  - 
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is  a  plane  wave  reflection  matrix  referenced  to  the  ground  and  represents 
reflection  from  everything  below  ground  level  with  vacuum  above  it  (i.e.,  the 
elements  of  R(0)  are  Fresnel  reflection  coefficients) .  Standard  notation  is 
used  in  the  sense  that  the  first  subscript  of  the  elements  of  R  and  R  denotes 
the  polarization  of  the  incident  wave  (  n  for  TM,  l  for  TE)  and  the  second 
subscript  denotes  the  polarization  of  the  reflected  wave.  Equation  (2)  repre¬ 
sents  two  equations.  Writing  out  the  lower  of  these  gives 


Solving  for  the 


ratio  ed  (o)/hd  (o) 
ym  ' '  ym v  ' 


from  che  above  yields 
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(6) 


Next  it  is  observed  that  the  total  modal  fields  at  the  ground  are  by  virtue  of 
the  definition  of  R 


h  (o)  =  (1  +  n R |, ( 0  ))hd  (o)  (7) 

yrii  '  '  1 '  nr  /  ym'-  ' 

e  (o)  =  (1  +  i R, (8  ))ed  (o)  (8) 

ym v  '  '  ~  x  N  m '  '  ym  '  • 

By  dividing  Equation  (8)  by  Equation  (7)  and  introducing  Equation  (6)  it  fol¬ 
lows  that 
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(9) 


evm(o)  (l  +  iRi(V))  n*±i0m)  „R,I  <*m) 
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The  factor  in  Equation  (9)  which  involves  products  of  the  elements  of  R  and  R 
is  a  true  global  parameter  (i.e.,  independent  of  source  and  receiver  geometry; 
and  is  independent  of  the  reference  height  chosen  for  those  elements.  This 
factor  will  be  referred  to  as  the  modal  polarization.  The  other  factor  in 
Equation  (9)  does  depend  upon  the  reference  height  chosen  for  R  and  R.  More 
generally 


e  ( z ) 
vm v 

h  (z^ 
ym 


1  +  iMVl 

1  +  "MVjz 


i.MV 

-  -  )|Rr(0 


where  the  subscript  z  denotes  the  reference  height  for  the  R  reflection  ele¬ 
ments.  Equation  (10)  gives  the  admixture  of  TE  and  TM  components anywhe re 
within  the  waveguide.  It  is  clear  that  the  modal  polarization  is  proportional 
to  the  conversion  coefficient,  (i ,  so  that  a  vertical  transmitter  which 
launches  only  TM  wave  will  excite  modes  of  hybrid  structure.  In  the  present 
case  interest  is  principally  in  the  effect  of  the  TE  admixture  at  the  ground. 
Since  at  the  ground,  LRL  =  -1  +  e  with  |e|«l,  the  factor  (1  +  j_Rx(0m))z=o 

constrains  the  field  to  be  very  small  relative  to  h^  at  the  ground. 

Nevertheless,  it  will  be  shown  in  Section  IV  that  the  longitudinal  component 
of  the  scaled  rf  magnetic  field, 


,  <9E 

u  =  -I  — Z 
x  ik  3Z 
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can  cause  significant  DOA  error  under  nighttime  propagation  conditions  when 
nRj_  is  the  largest.  This  is  particularly  true  in  null  regions  of 


H  .  Included  in  the  results  of  Section  IV  are  the  mode  sums 
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(13) 


evaluated  at  the  ground  for  ground  based  vertical  transmitters.  It  is 

remarked  that  h  and  e  are  simply  linear  combinations  of  Airy  functions  (or 
ym  ym  r  •'  J 

linearly  related  functions)  and  so  the  derivative  term  in  Equation  (13)  is 

readily  calculable.  Formulas  for  the  height  gains  are  given,  for  example,  in 

reference  6. 

If  the  mode  sums  are  written  as 

Hx  “  a^cos(wt+5^)  (14) 

Hy  -  a2Cos (wt+52 ) ,  (15) 

the  end  point  of  the  magnetic  intensity  traces  out  an  ellipse  (see  Figure  1) , 
the  properties  of  which  are  immediately  obtainable  from  reference  9 . 


ng.  1  MAGNETIC  INTENSITY  ELLIPSE. 
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In  particular  the  angle  ip  in  Figure  1  is  given  by 


4>  -  1/2  tan  \tan2acos5)  (15) 


where 


tana 


o  <  a 


7T 

2’ 


(16) 


Recalling  that  x  -  z  is  the  plane  of  propagation  with  x  the  range  coordinate 
and  that  y  is  the  transverse  coordinate  in  a  right  handed  coordinate  system 
the  DOA  error  is 


X  “  */2  -  ■  (17) 

Positive  x  implies  arrival  from  the  second  quadrant  of  Figure  1,  whereas 
negative  x  implies  arrival  from  the  third  quadrant.  In  section  IV  examples  of 
mode  sum  plots  for  and  H.  are  given  along  with  plots  showing  the  DOA  error 

y  x 

for  paths  examined  in  reference  2. 


TIT.  BRIEF  PROGRAM  DESCRIPTION 


GENERAL  COMMENTS 


The  program  discussed  in  this  section  is  a  modification  of  the  program  dis¬ 
cussed  In  reference  3  and  can  be  used  to  treat  either  horizontally  homogeneous 
or  inhomogeneous  guides.  To  handle  horizontal  inhomogeneities  the  waveguide 
is  divided  into  a  series  of  slabs  as  shown  in  Figure  2  below: 


Figure  2.  Mode  conversion  model. 


In  particular,  Figure  2  shows  the  convention  for  slab  numbering  (which  is 
reversed  from  that  of  reference  3)  and  a  possible  configuration  for  a  ter¬ 
minator  type  lateral  inhomogeneitv.  The  transmitter  is  at  the  origin  of  the 
coordinate  system  and  the  value  x-,  ,  which  is  not  shown,  is  not  used  in  the 

calculations  but  must  be  assigned  a  value  in  namelist  DATUM  because  of  array 
dimensioning.  Two  options  are  available  with  the  program.  For  each  slab  and 
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for  each  mode  the  ground  eigenangles  must  be  supplied  along  with  four  inde¬ 
pendent  quantities  from  which  electric  point  dipole  excitation  factors  are 
determined.  These  quantities  are  supplied  by  the  waveguide  program 
"MODESRCH" .  Unlike  reference  3,  the  number  of  modes  can  vary  from  slab  to 
slab  in  the  present  program.  Also,  the  quantity  F0FR  defined  by  Equation  (5) 
of  reference  3  (also,  equivalent  to  Equation  (9)  of  the  present  study)  is  cal¬ 
culated  internally  in  the  present  program. 

Mode  sum  plots  for  H  ,  Hx  and  the  DOA  error  can  be  generated  as  a  function  of 

range  x  for  a  fixed  position  of  the  lateral  inhomogeneity  or  they  may  be 
plotted  as  a  function  of  the  tranmitter  terminator  distances,  D,  for  a  fixed 
receiver  site.  The  latter  option  can  be  used  only  when  the  ground  electrical 
properties  and  geomagnetic  field  magnitude  and  orientation  can  be  taken  as 
constant  over  the  path.  The  plots  can  be  generated  at  an  arbitrary  height 
within  the  waveguide  for  electric  dipole  sources  of  arbitrary  orientation  and 
height  within  the  waveguide  (see  Figure  3) . 


z 


Fijure  3.  Dip  de  M  orientation  within  the  waveguide  where  7 
is  the  inclination  ar.J  ?  the  azimuthal  orientation. 
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B.  DESCRIPTION  OF  INPUT 


All  input  to  the  mode  conversion  program  is  given  in  a  data  file  on  the  stan¬ 
dard  input  unit.  A  listing  of  sample  input,  shuwin6  the  data  file  setup  for  a 
range  calculation  for  a  horizontally  homogeneous  guide  is  given  on  page  12. 

An  input  file  setup  for  a  moving  terminator  is  given  on  pages  13  through  15. 

A  sample  input  for  a  range  variation  in  a  laterally  inhomogeneous  guide  would 
be  similar  to  the  sample  given  in  reference  3. 

There  are  two  parts  to  the  input.  The  first  part  is  read  in  by  means  of  a 
Fortran  77  namelist- directed  read  statement  with  a  group  name  of  DATUM.  The 
following  variables  and  arrays  should  be  specified  in  the  NAMELIST  input. 

RHOMAX  -  maximum  range  in  km  for  which  calculations  are  performed. 

RHOMIN  -  minimum  range  in  km  for  which  calculations  are  performed. 

DELRHO  -  horizontal  increment  in  km  for  which  successive  field  strengths  are 
computed . 

NRSLAB  -  Number  of  slabs  in  the  model. 

NTMAX  -  Number  of  times  the  tiansmi tter  terminator  separation  is  incremented. 

DELTAX  -  Distance  in  km  by  which  transmitter  -  terminator  separation  is  incre¬ 
mented  . 


11 


SAMPLE  INPUT  FOR  LATERALLY  HOMOGENEOUS  GUIDE 

$ da turn 

rhomin— 25. ,  rhomax-10000 . .delrho— 25. , 
deltax-0 . , ntmax-1 , 
nrslab-1 , 
xval-0 . , 

xmin-0. , xmax-10000 . ,xtic=1000. , 

ymin--90 . ,ymax=90 . ,ytic=10 . , sizex=5 . , sizey=6 . , 

gamma-45 . , phi-45 . , 

zr-10 . , zt=10 . , 

iprnta-0 , 

ipltop=2 , 

iplflg-2, 

if irst-1 , 

last-1 , 

$end 

Lualualei  to  San  Diego  beta- . 7 , hprime-87 . 0 

R  .000  F  23.4000  A  59.000  C  40.000  M  .413E-04  S  4.000E+00  E  81.0  T  87.0 

1  89.97682  -5.920991  1 . 17507458E-04-1 .42930751E-03-2 . 91159674E-12  1 . 56223735E-13 

2  89.97682  -5.920991  5 . 15993896E-08  2 . 86593256E- 08  1 . 00046158E+00  6 . 63018107E-01 

1  89.82555  -5.574832  2 . 17080335E-04- 1 . 72104372E-03 - 5 . 10516109E- 12 - 1 . 08497217E- 12 

2  89.82555  -5 . 574832-6 . 85928470E-08-5 . 32394253E-08  1 . 00075912E+00  6 . 63700819E-01 

1  87.88632  - . 224351-2 . 35030573E-04-2 . 65855566E-02 - 1 . 17070893E- 12- 1 . 08537094E- 12 

2  87.88632  -.224351  1 . 22134551E-07  1 . 42263872E-07  1 . 00533319E+00  6 . 68529391E-01 

1  84.87724  -.264002  9 . 20612714E-04- 1 . 01905130E-03 -9 . 30642854E- 11- 3 . 59403197E-12 

2  84.87724  -. 264002-2 . 07534725E-07 -2 . 50064232E-07  1 . 00777233E+00  6 . 71742201E-01 

1  81.56534  - .157771-9.38259356E-04-2.64533665E-02-7.21719646E-12-7.26697782E-12 

2  81.56534  -.157771  3.05091049E  07  3 . 59546164E-07  1 . 01319230E+00  6 . 78170621E-01 

1  79.29093  -.254342  1 . 24309910E-03- 1 . 29618787E-03 -2 . 71744988E-10-1 . 65220084E- 12 

2  79.29093  -  .  254342 -4 . 04487480E-07 -4 . 83902454E-07  1 . 01797378E+00  6 . 84924006F.-01 

1  76.56423  - . 180111-1 . 40988897E-03-2 . 33948622E-02 -2 . 79152760E- 11 -  2 . 70179695E-11 

2  76.56423  -.180111  5 . 70840825E-07  6 . 39522682E-07  1 . 02576578E+00  6 . 94986641E-01 

1  74.45299  -.317182  1 . 86981435E-03 - 1 . 98888732E-03-5 . 12932530E-10  1 . 35284353E- 11 

2  74.45299  - . 317182 -6 . 98896315E-07- 7 . 95034907E-07  1 . 03191864E+00  7 . 05408871E-01 

1  71.81008  - .208421-2.10860302E-03-2.16372982E-02-7.77475515E-11-7.09266951E-11 

2  71.81008  -.208421  9 . 32243722E-07  9 . 67916662E-07  1 . 04210842E+00  7 . 20323980E-01 

1  69.73547  -.405862  2 . 68813735E-03 -2 . 99302535E-03- 8 . 04118716E-10  5 . 29969506E- 11 

2  69.73547  -. 405862 - 1 . 09387440E-06- 1 . 15510568E-06  1 . 04884636E+00  7 . 35498965E-01 

1  67.09178  - .237161-2.93348846E-03-2.01747213E-02-1.77707127E-10-1.50172735E-10 

2  67.09178  -.237161  1 . 39038355E-06  1 . 30640410E-06  1 . 06097567E+00  7 . 56687820E-01 

1  64.97400  -.520412  3 . 63010028E-03-4 . 30195918E-03 - 1 . 13163612E-09  1 . 28242014E- 10 

2  64.97400  -. 520412-1 . 59059425E-06-1 . 52464986E-06  1 . 06673014E+00  7 . 78390706E-01 

1  62.30975  - .265031-3.78178176E-03-1.87656302E-02-3.51571799E-10-2.72375650E-10 

2  62.30975  -.265031  1 . 93094775E-06  1 . 60819002E-06  1 . 07958007E+00  8 . 07601869E-01 

y 

y 
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SAMPLE  INPUT  FOR  MODELING  A  MOVING  TERMINATOR 

$ da turn 

rhomin-3821. , rhomax-4166  . , delrho-345 . , 
deltax-25. ,ntmax=261, 
nr slab-9 , 

xval— 1525. ,-1500. ,-1312.5,-1187.5,-1062.5,-937.5,-812.5,-687.5,-500.  , 

xmin--1500. , xmax-5000. , xtic-500. , 

ymin-- 30  .  ,ymax-50 .  ,ytic-10 .  ,  sizex-7 . 82  ,  sizey-4 . 92  , 

gamma-45, , phi=270. 

zr-15 . , zt=15 . , 

iprnta-0 , 

ipltop-1 , 

iplflg-1, 

if irst-1 , 

last-1 , 

$end 

Hprime-86  to  Hprime=70 

R  .000  F  21.7940  A  58.500  C  39.000  M  .431E-04  S  4.640E+00  E  81.0  T  86.0 

1  89.95893  -5.614791  8 . 90918891E-05-2 . 56474130E-03-3 . 66113606E- 12 -1 . 41053103E-13 

2  89.95893  -5.614791  7 . 23541476E-08  5 . 82935336E-08  1 . 04440999E+00  3 . 08963150E-01 

1  89.73230  -5.098552  8 . 13409977E-04- 3 . 18036205E-03-7 . 20839404E-12-1 . 23412988E-12 

2  89.73230  -5 . 098552-9 . 77022978E-08 - 1 . 11584257E-07  1 . 04662442E+00  3 . 09465617E-01 

1  86.46790  - . 276901-6. 97296404E-04-2. 96539459E-02-2 . 31411700E-12-2 .42465618E-12 

2  86.46790  -.276901  1 . 52632509E-07  2 . 58247184E-07  1 . 06115580E+00  3 . 07544500E-01 

1  83.71133  -.318972  1 . 78535702E-03- 1 . 91355264E-03-1 . 02256321E-10-2 .46968201E-12 

2  83.71133  - .318972-2. 50906027E-07-4. 20881236E-07  1 . 07172537E+00  3 . 07309061E-01 

1  80.24526  - .275711-2.00039987E-03-2.66215149E-02-1.59478784E-11-1.70393977E-11 

2  80.24526  -.275711  3 . 87197304E-07  6 . 31316766E-07  1 . 09389603E+00  3 . 06244135E-01 

R  .000  F  21.7940  A  58.500  C  39.000  M  .431E-04  S  4.640E+00  E  81.0  T  84.0 

1  89.93198  -5.353161  1 . 86523030E-04- 3 . 80378054E-03- 3 . 72650522E- 12-2 . 66069854E-13 

2  89.93198  -5.353161  9 . 06502038E-08  7 . 63756844E-08  9 . 70103979E-01  2 . 90322542E-01 

1  89.67323  -4.842402  9 . 74733557E-04- 3 . 26375384E-03-9 . 61524870E-12 -1 . 76274479E-12 

2  89.67323  -4 . 842402 - 1 . 15998247E-07 - 1 . 39092066E-07  9 . 72432077E-01  2 . 91206360E-01 

1  85.93787  - . 315181-6 . 77048694E-04- 3 . 09771523E-02-2 . 07085325E-12-2 . 71785220E-12 

2  85.93787  -.315181  1 . 47816920E-07  2 . 84299375E-07  9 . 89267349E-01  2 . 90921926E-01 

1  83.22681  -.388222  1 . 78511837E-03-1 . 55677041E-03-1 . 17126780E- 10-4 . 95431560E-12 

2  83.22681  -. 388222-2 . 37927182E-07-4. 57677174E-07  1 . 00180018E+00  2 . 92327493E-01 

1  79.73028  - .330761-2.05988088E-03-2.73547117E-02-1.47570515E-11-1.97417777E-11 

2  79.73028  -.330761  3 . 79882323E-07  6 . 98669737E-07  1 . 02771771E+00  2 . 94379383E-01 

R  .000  F  21.7940  A  58.500  C  39.000  M  .431E-04  S  4.640E+00  E  81.0  T  82.0 

1  89.88525  -5.073801  4 . 76708548E-04-5 . 71369240E-03-3 . 24168990E- 12 -3 . 47774489E- 13 

2  89.88525  -5.073801  1 . 07369090E-07  9 . 53978301E-08  8 . 63476753E-01  2 . 74763733E-01 

1  89.59141  -4.554722  9 . 79591976E-04- 2 . 88872421E-03- 1 . 31660897E-11-2 . 77030027E-12 

2  89.59141  -4. 554722-1. 32362615E-07-1 . 66220417E-07  8 . 66081953E-01  2 . 75995642E-01 

1  85.41425  - .371091- 5. 08478784E- 04-3. 22124176E -02-1. 73904988E- 12-2. 72687819E-12 

2  85.41425  -.371091  1 . 43651306E-07  3 . 02742137E-07  8 . 85840416E-01  2 . 76717037E-01 

1  82.72205  -.487472  1 . 59731111E-03- 1 . 18575653E-03- 1 . 34468672E- 10-9 . 75490175E- 12 

2  82.72205  - .487472-2 . 24028838E-07-4. 83897225E-07  9 . 00848508E-01  2 . 79662967E-01 

1  79.18263  - .403351-1. 88347534E-03-2 . 82151829E-02-1 . 27613467E-11-2 .08172975E-11 

2  79.18263  -.403351  3 . 71068296E-07  7 . 55888038E-07  9 . 31417108E-01  2 . 84143299E-01 

R  .000  F  21.7940  A  58.500  C  39.000  M  .431E-04  S  4.640E+00  E  81.0  T  80.0 

1  89.82103  -4.781621  1 . 05937861E-03 - 8 . 39334913E-03- 2 . 11253519E- 12 -2 . 83353229E- 13 

2  89.82103  -4.781621  1 . 13315721E-07  1 . 01717518E-07  7 . 26135790E-01  2 . 78738737E-01 

1  89.46140  -4.205672  7 . 41320197E-04- 1 . 97134237E-03 - 1 . 82106691E- 11-4 . 62774090E- 12 

2  89.46140  -4. 205672-1. 38748817E-07-1 . 77043150E-07  7 . 30000734E-01  2 . 80761808E-01 

1  84.88551  - . 447061-1. 71756066E-04- 3 . 33595462E-02 - 1 . 36530576E- 12  -  2 . 37677196E- 12 


2  84.88551  -.447061  1 . 44529494E-07  3 . 05609575E-07  7 .49675095E-01  2 . 79543728E-01 

1  82.17252  -.618432  1 . 24706491E-03 - 7 . 74172600E-04-1 . 54974325E- 10- 1 . 78294827E-11 

2  82.17252  -. 618432-1 . 94097282E-07-4 . 85529313E-07  7 . 95087934E-01  2 . 63119757E-01 

1  78.59708  -. 499461- 1.50697830E-03-2. 91476566E-02 -9. 99454252E- 12-2. 02260205E- 11 

2  78.59708  -.499461  3 . 56034121E-07  8 . 13086274E-07  7 . 94667065E-01  2 . 58789837E-01 

R  .000  F  21.7940  A  58.500  C  39.000  M  .431E-04  S  4.640E+00  E  81.0  T  78.0 

1  89.74937  -4.485971  1 . 78524549E-03 - 1 . 14039388E-02 -9 . 32397931E- 13- 1 . 85739824E- 13 

2  89.74937  -4.485971  9 . 78692469E-08  8 . 47686010E-08  5 . 76031327E-01  3 . 11028749E-01 

1  89.24081  -3.753632  4 . 48389386E-04-9 . 48478992E-04-2 . 46174071E- 11-7 . 50462800E- 12 

2  89.24081  -3. 753632-1. 24890576E-07-1 . 59517882E-07  5 . 78858852E-01  3 . 12637269E-01 

1  84.34246  -.546751  2 . 11181687E-04- 3 . 44251283E-02- 8 . 89651688E- 13 -1 . 74506330E- 12 

2  84.34246  -.546751  1 . 43204801E-07  2 . 80839743E-07  6 . 02452695E-01  3 . 11667144E-01 

1  81.56769  -.801902  8 . 28500022E-04- 5 . 05273812E-04- 1 . 80112383E- 10- 3 . 04999290E- 11 

2  81.56769  -. 801902-2 . 07465988E-07 -4 . 59239033E-07  6 . 22810721E-01  3 . 15373600E-01 

1  77.96826  - . 635241-7 . 57292029E-04- 3 . 02124228E-02-7 . 16761546E- 12  - 1 . 57295739E- 11 

2  77.96826  -.635241  3 . 67441686E-07  7 . 59492764E-07  6 . 61270738E-01  3 . 18266243E-01 

R  .000  F  21.7940  A  58.500  G  39.000  M  .431E-04  S  4 . 640E+00  E  81.0  T  76.0 

1  89.67842  -4.195361  2 . 45779543E-03 - 1 . 41851194E-02 - 3 . 12425080E-13- 1 . 36429022E- 13 

2  89.67842  -4.195361  7 . 04787269E-08  5 . 99233658E-08  4 . 21555966E-01  3 . 88745427E-01 

1  88.90660  -3.171932  2 . 70018878E-04-3 . 39354476E-04- 3 . 21954893E- 11-1 . 11604580E- 11 

2  88.90660  - 3 . 171932 - 9 . 63160147E-08 - 1 . 28835026E-07  4. 29769784E-01  3 . 76125246E-01 

1  83.78646  -.656921  5 . 78613428E-04- 3 . 53961140E-02-4 . 17128761E-13- 1 . 14611727E- 12 

2  83.78646  -.656921  1 . 34711897E-07  2 . 33643163E-07  4 . 48909104E-01  3 . 88394594E-01 

1  80.88498  -.976642  4 . 98450419E-04-2 . 34734369E-04-2 . 11642412E- 10-4 . 63897462E- 11 

2  80.88498  -. 976642 - 1 . 95438730E-07 - 3 . 96136898E-07  4 . 70917970E-01  3 . 90501976E-01 

1  77.29854  -. 793811-2 . 86844370E-05-3 . 13108861E-02- 3 . 75179878E- 12- 1 . 10646362E- 11 

2  77.29854  -.793811  3 . 55859243E-07  6 . 66429230E-07  5 . 09032011E-01  3 . 89449626E-01 

R  .000  F  21.7940  A  58.500  C  39.000  M  .431E-04  S  4.640E+00  E  81.0  T  74.0 

1  89.60178  -3.902511  3 . 09083378E-03-1 . 67904254E-02 -9 . 43353129E-14-8 . 49416295E-14 

2  89.60178  -3.902511  4 . 56988438E-08  3 . 96695583E-08  2 . 91040480E-01  5 . 15212834E-01 

1  88.40370  -2.497722  1 . 53187531E-04- 1 . 08900967E-04-4 . 13913175E- 11- 1 . 51935877E- 11 

2  88.40370  -2 . 497722 - 7 . 22185121E-08-9 . 34580413E-08  2 . 95880765E-01  5 . 15322745E-01 

1  83.21597  -.773371  8 . 99333565E-04-3 . 62679623E-02 -1 . 01535987E- 13- 7 . 34521602E- 13 

2  83.21597  -.773371  1 . 13296082E-07  1 . 78915556E-07  3 . 16028774E-01  5 . 09744406E-01 

1  80.15169  -1.113712  2 . 94007245E-04-7 . 68069804E-05 -2 . 50641397E- 10-6 . 35519415E- 11 

2  80.15169  -1. 113712-1. 70118525E-07-3 . 16608464E-07  3 . 39071542E-01  5 . 05254328E-01 

1  76.58773  -.971591  6 . 23005268E-04- 3 . 23789045E-02 - 1 . 09145711E- 12- 7 . 29663465E- 12 

2  76.53773  -.971591  3 . 13271016E-07  5 . 29444264E-07  3 . 73405904E-01  5 . 06903291E-01 

R  .000  F  21.7940  A  58.500  C  39.000  M  .431E-04  S  4 . 640E+00  E  81.0  T  72.0 

1  89.50670  -3.593371  3 . 73883406E-03- 1 . 93957984E-02 - 2 . 80016150E- 14-4 . 41833222E-14 

2  89.50670  -3.593371  2 . 83041288E-08  2 . 61669086E-08  2 . 04742640E-01  6 . 65039361E-01 

1  87.66309  -1.927862  8 . 04060764E-05 - 3 . 38098434E-05 - 5 . 22675826E- 11- 1 . 98142700E- 11 

2  87.66309  -1 . 927862 - 5 . 04189011E-08 - 6 . 68672868E-08  2 . 10287109E-01  6 . 64982557E-01 

1  82.62367  -.904771  1 . 22795266E-03- 3 . 70658040E-02  2 . 93056933E- 14-4 . 71066925E- 13 

2  82.62367  -.904771  8 . 56292814E-08  1 . 32001347E-07  2 . 27945834E-01  6 . 60904884E-01 

1  79.40465  -1.226292  1 . 72982021E-04- 1 . 48205954E-05 -2 . 95881680E- 10- 8 . 13693696E- 11 

2  79.40465  - 1 . 226292 - 1 . 35712256E-07-2 . 37786438E-07  2 . 47396663E-01  6 . 59889698E-01 

1  75.83375  -1.174621  1 . 18812581E-03-3 . 33632752E-02  2 . 17520067E- 13 -4 . 73121412E- 12 

2  75.83375  -1.174621  2 .44619088E-07  4 . 01201845E-07  2 . 78375298E-01  6 . 56392694E-01 

R  .000  F  21.7940  A  58.500  C  39.000  M  .431E-04  S  4.640E+00  E  81.0  T  70.0 

1  89.37548  -3.251751  4.42431914E-03-2 . 21310500E-02-8 . 24005764E- 15 - 1 . 95039839E-14 

2  89.37548  -3.251751  1 . 75193087E-08  1 . 72535923E- 08  1 . 63431019E-01  8 . 03890228E-01 
1  86.81053  -1.622792  3 . 89222732E-05 - 1 . 11819172E-05 - 6 . 52232504E-  11  -  2 . 57851605E- 11 
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c 


2  86.81053  -1. 622792-3. 47190507E-08 -4 . 73781974E-08  1 . 69330418E-01  8 . 04408669E-01 

1  81.99325  -1.059281  1 . 60216109E-03- 3 . 78431045E-02  4 . 95666469E-14- 2 . 86521430E-13 

2  81.99325  -1.059281  6 . 19566336E-08  9 . 75383330E-08  1 . 85555279E-01  8 . 03614676E-01 

1  78.64407  -1.357022  9 . 65972431E-05-4 . 17938452E-07 -3 . 47162410E- 10- 1 . 02328881E- 10 

2  78.64407  - 1 . 357022 - 1 . 01884517E-07 - 1 . 78062535E-07  2 . 05252454E-01  8 . 05077076E-01 

1  75.02145  -1.414181  1 . 72621966E-03 - 3 . 42926793E-02  5 . 21807912E-13-2 . 98834373E-12 

2  75.02145  -1.414181  1 . 78861782E-07  3 .04132669E-07  2 . 34394595E-01  8 . 06276321E-01 
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XVAL 


IFIRST 


LAST  - 


IPLTOP 


XMIN  - 


XMAX  - 


XTIC  - 


YMIN  - 


Horizontal  position  in  km  of  boundaries  between  adjacent  slabs.  These 
are  denoted  by  the  x^'s  in  Figure  2.  Note  that  XVAL  can  be  negative 
and  that  it  is  dimensioned  for  25. 

-  Is  set  to  1  in  the  first  set  of  NAMELIST  input.  If  more  than  one  set 
of  input  is  used  set  IFIRST  =  0  in  the  second  set. 

Is  set  to  1  in  the  last  set  of  NAMELIST  input.  If  the  user  has  re¬ 
quested  plots  this  causes  the  end  of  file  to  be  written  on  the  plot 
tape  . 


-  Plotting  option  flag.  If  IPLTOP  =  1,  mode  sums  for  the  scaled  rf 
magnetic  field  components  H^_,  and  the  DOA  error  are  plotted  as  a 
function  of  transmitter- terminator  distance  at  a  given  field  point. 


If  IPLTOP  =2,  and  the  DOA  error  are  plotted  as  a  function  of 

range . 


Minimum  value  of  x  on  plot  (x  either  range  or  transmitter- termi nator 
distance) . 


Maximum  value  of  x  on  plot  (x  either-  range  or  transmitter- terminator 
distance )  . 

Tic  mark  interval  on  x  axis. 


Minimum  value  of  v  on  plot  (y  is  H^,  and  DOA  error).  Normal  setting 
is  -90. 
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YTIC  -  Tic  mark  interval  on  y  axis. 


SIZEX  -  Size  of  x  axis  in  inches. 


SIZEY  -  Size  of  y  axis  in  inches. 


GAMMA  -  Dipole  orientation  relative  to  z  axis  (see  Figure  3) . 


PHI  -  Dipole  orientation  relative  to  x  axis  (see  Figure  3) . 


ZT  -  Transmitter  altitude  in  km. 


ZR  -  Receiver  altitude  in  km. 


INTFLG  -  Printing  option  flag.  INTFLG  must  be  set  to  1  if  printout  of  height 
gain  integrals  is  required. 

IPRNTA  -  Printing  option  flag.  IPRNTA  must  be  set  to  1  if  printout  of 
generalized  mode  conversion  coefficients  is  required. 

IPLFLG  -  if  IPLFLG  f  0  plots  are  determined  by  IPLTOP.  If  IPLFLG  -  0  no  plots 
are  generated. 
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The  v's  or  n's  appearing  at  the  end  of  the  input  files  on  pages  12  and  15  are 
in  response  to  the  questions  appearing  at  the  end  of  the  output  files  on  pages 
27  and  37. 

The  second  part  of  the  input  contains  NPUNCH  =  1  output  from  "MODESRCH"  .  The 
output  file  contains  the  frequency  in  kHz  and  the  following  slab  data. 

SIGMA  -  Ground  conductivity  in  Si/m.  SIGMA  is  dimensioned  for  25  slabs. 

EPSR  -  Relative  permittivity  of  the  ground.  EPSR  is  dimensioned  for  25  slabs. 

TOPHT  -  A  control  height  in  km  for  the  upper  limit  on  the  height  gain  in¬ 
tegrals  (see  reference  6).  TOPHT  is  dimensioned  for  25  slabs. 

IDPLOT  -  Literal  constant  tip  to  68  characters  which  is  printed  on  plots 
produced.  For  example,  in  the  NAMELIST  input  data,  place  the  card 
IDPLOT  =  "CUTLER  TO  SAN  DIEGO  PATH". 

The  NPUNCH  =  1  output  file  from  "MODESRCH"  also  contains  for  each  mode  the 
ground  eigenangles  along  with  four  auxiliary  quantities  (called  Tl,  T2 ,  T3  and 
T4 .  see  reference  3)  for  calculating  excitation  factors. 

C.  OUTPUT 

The  sample  output  shown  on  pages  1 C)  through  27  corresponds  to  the  input  on 
page  12.  Figure  4  is  the  plot  output.  The  principal  hardcopy  output  begins 
with  RHO-KM.  MAG ( HY  i  - DL . HAG. HX ) - DB  and  DOA  ERROR -DEG.  RHO  is  the  range 
MAC ( HY )  | H_  | .  MAG (MX)  -  | H  |  in  DB  above  Ipv/m  and  DOA  is  the  direction  of 
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rn^fNCN<NOvo<Nr3'^'iocr\oJoricg'<rU)r'<NrricriLr)ror^r'CT'(Nir'OL?)r^Cor~<-Htr)in>-('>a,'^rriO(Nji-i1HO[^(Nr'CTicr>c\ju)r' 

vOCVD)OO^vOCva>QfNr^vO^TriDVOvOO^OOC'\i'*},Hr\)r','><7>vOCV(T>r-i(7ir''fNJVOr^COCOr'">rO'34CO('\jr\lrHroovvOr\)C\JOcovOCO'3* 

TriD^OiDav«H<HOr^rocorHc\iDirHcocN^cNiiDricoHr^vocNHawiDiDrococNiDr~cor~iDC'viav'"rcoi--iroiDr^ovo,sr'3,voo 

id  vo  d- oo  oo  oo  cv  ci  crv  cc  co  t'' vo  id <n  rH  i  01  ^  cp  vo  o  aw  HeHcvw  id  vor~  rococo  coco  coco  cot-'t^vo  void  ^nn  cm  cndj  cm  o 

I  I  I  rj-IDJ  (  |  | 


ooQv<N^vor~cjir^rH>-ifMCMovoiDvoco'i>(Tivooo^Hdr^O'^rt^o^}‘CNooc\coiDiD'^,^j,vovDc'javc3vvonvooor'L 
■'TnocNvpcrKN^'^T'rHcocDO^rDr'avoo'fiDrHCNt-iocvicocMiDcnrir'tDTy^'inr^c^oooocotHiDi-iavooavncoS'^TM 
cm  00  ea  O  t-- o  oo  m  o  a\  n  ro  n  cm  cv  o  <d  ^  eg  ^  f'1  ^  ^  c"1  co  £t  vo  id  u"w  vo  vo  o - -  - -  ■  ~  ■ 


_  _  j  r- id  r-w  n  vo 

-  _  J  00  C0  iH  ID  i-H  O'  O  O  O'  CD  CD  t''  CM 

.  - ,  _.....  — - ,  .  - - , - >iDcv)<-iaico<-Havc'wioopcor''rwr''-CT>r-w 

r^coiDr'r'CN)TTiHcriOV'HoviDfnrncv)rHvpvovovovor^coor-iQrHiDi-iavcovorvivovoriiDi--irnr-(ro(^vOt^iDiHiDcoOf-ii-irHcTi 

cvjavvorgco^awcocoaw<HcvcococvO(-'W'OcoO(NLO[^onvoocor'-HiDco<-wvocoavooocrico[^-vO'*rrvir-ic'tv'LDc\i 

rMooo<j\avtooor'r'VovovoiDiDiDiDvovovovovo[v-r''r~r'OOQococna>cooddr4iHr-3(HrHDj(NDj<4rH<HrHrHrHiHddod 

^^^^ocoo(o<ororo(oco<ococicocn(oraroroororr)c')<'")C')r">r'">cor'‘)'‘T'3W'’3W,'3W''i‘'3W''3W'>aw''aw'sjw«qr^w*3> 
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r'OcMiDr'OcMiDr'0<MiD[^0(MiDr-orsuDr~oojiDt^ODHDr^o<MiDr'OCNiD[^ocMiDr'ODJiDr'OcNaDr^OfMiDr'0 

OJ'^C?f^r25T,3"3.‘'5Cir>lC?ir)LCl'£l^^:)vi:)r^r'r'r^(:9(5I5c2C^<^c^^oooO'-,'-<r-|'-|DiDjDjcMnmcoM'T'T^-<TlDiDiDtDvo 

idididididiDiDiDiDididididididididididididididldididldididvovovovovovovovovovovovovovovovovovovovovovovovovo 
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6625.00  40.06582  3.39362  .70249 


O'TiojMinr^^inroiQroinvo^'t'io^voortr'CNHinr^^oiriCTityi'SPrrLnincyifor'Vo^owJinr-iroocofn'a-cofNilStTivDH^vZ) 

m.no^cM.Hg}rHvoo>oc2cnawn'^'^u)Hr^^inm<nmvou?'£mo'«3,t^mr>CT\<T*cy'C'~Tr<T<rH®'00co(T>c\cov5ovo*j‘ 

Lncr(<Nrririf\ioocncn(NCncrirvi^rTfojcor)oo<Nr~fHvocNcO''rO|25<Nr'tHTi<r'Q&r'y3ir)^r<NrHcoLnr>jr^cvj^)(Tii-tcNior' 
r-r^cooococooocor^r^ooooco<riCTicr»cr>cocor'r'VOvoimn'f'a‘^,nncMrgrHooc5«-toon^riovoCDr'OococT»CTv&i00o<j» 


I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  -HrHrH  | 


cq^Q^(^rov£)q;r^(No0c5cgrnr-H^o^roLnoS'^,coricosOtNOi-(a>cr'r^O'>toriCT\(viri(NrH'tj'0(Ti'jD^Hcor)(TiLOcoiri 

^t>i)Oor'fNojinr'r'OocMg^noor;Lr)com^oO'T3p-ic^Orv<fnr^tnvOrHc5r^vonvDLf)rH^Dv£>TrHr>)MCMmr)ir)a<o 

t^,<H^cor^-cocoir)^rrocorr)OLn^,CTiQrHrofNcrir'OocQror-iTHCTiv£)<rior'rHCQinvDcQOto<bror^<rio>r4oo’a'r>rivDoor'ai 

rHrHcr>rocr>ogOTU5'^'CT<o>jcif^ojr^o>i/it^y3r^r-icofNj(5irii/,)a)<Nri’trr~<viovicOfV(^a5’g,<viv5f\>ooir)^'voeorvir'^ffn'i,oor' 

M)r~^>tOi-i'i5c3TrcocNOOi-i'^,vocooo(nr'iDi-(<\jvoOfoinr^r^LnHvoor~^D''i,<MOtv'ncy»nr^o<Mn'j’OfMOr^no' 


ononriNNH 


I  t  pH  rH  cH  *— I  rH  rH  rH  t 
I  I  I  I  I  I  I 


pHcNM'a,rfimr)«)v£>r^r'r'r^r>‘r^r'VOvoio 


^■niHcnn<Tir^cocO”3,crirHfvjorHOr'Ocgrvi^or'C7iTf'i-tTfiHrrunvooir)r~co'£>(7icnvOr-ityiOTj*voir)rM(McrunrN^r«3*(Tvo' 

<j\r^^CTvpHfo-^‘r,i'>3,(N<Nrorir>ir)vovO'j,kDHrnr^<j\'iir-inri«a'^r'^)iH{NO^rco^f(NC~-tNinvoinrHOr'iOCNrHrHco'>3,coCT\ 

cot^ovDrrrHM5^3,pHinor^ONOrH^rcc^£)ojH(Ncomx)pHnr''i)'X)^^3(NncoCorHCMcoQn^},iTir^oo<M<-tvOf-tr-ioiO'-ir~ 

drJl°^r^CT'.c:cl[v'CrCT't^C^<:Dvor1$<^'7,Cl!£)coa:i^)Of^'H®'!:i‘0>^fOMnrHr"THr^iHi7>vouSrHin3if)rH^fvovoi?ir)CN'-' 

coincN(j>uii-ir~friOMr)CNiococoo\or\jiocooojTryD(X)(TiOo<HpH(N<NnTj'Tj<in\>ovor~c^t^t^r~voinTrfHHCT>r'if)nr4o\ 

<Ticncr\cocooor^t^vcvovovoiritnir)vou50^r^r^[^t^r^r'-cococococoeocococooococococococococococooocDt''r'-C''-r~r'-'0 

rHfotorirnMrHfH^rooMcHo^ricoornronfnoonroricnrHfnfHcncnrHririrHroofnrHfOPHfnrHrorirHririoriri 


88888888888888888888888888888888888888888888888888888 
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r;  gj  r- oj  cri  oo  in  vo  n  n  h  o  H  co  co  Q  Q  ir>  vo  rr-i  <N  o  vo  <n  rH  tt  cn  <h  r^vo  n  irun  cm  om>- rH  o  oq  vd n  co  rH 

inoc\igwcT>vooc^n'>jvoro^t^in>HOixior-'<HHcomc^(N>-4^ot^vO'r'^om^mQOcNCVcHmrvjiHmt^®vS'OOr''in 
co>H30vOin«H-<jr}rHinrvioooH<\i^cvcriiHr^in'*r'srrNOinc\)(N!vocNcrivo<ot^co'?i7icriir><NrMnr~ncNvDrHor~ovocri 
cvi  r' o  co  in  co  cvi  vo  ^-i  co  o  rH  oo  in  cvi  gv  i^-  >o  ^  co  <n  <-(  q  <-i  f\j  g- vp  co  r-i  ro  vo  cti  <-i  co  in  vo  vo  vo  tn  n  o  vo  <n  i-i  o  in  o  in 

iJvcOCOr^vOinm^rfcOCOCOOlCMCNiHHHOOOOOOOOOOOOlDi-lHrHrHlMCMCsJlNCMrvJCNJCNJfvjr-lr-loOOi-IHOOOJCN 

r  r  r  r  r  r  r  r  r  r  r  •*  r  r  r  i*  r  i*  r  i*  i*  r  i*  r  r  i . r  r  r  r  r  i* 


cocovoinocsiiHc»<HovoincoiHr^vocor^cocor^coir»r'avovr-i(Nr-i®ocMir)r'i-Hvo(rivovoiNOOvo'3,coaM,ococor'Voir)rHin 

incocoav^,coor'0'^‘0'3,cr>i^vocococvj^rvovo'jr^^i,riir^'3,coinoojovoi-(r'coor^cM^},0'd‘i-ivOcooMX>’Tcvjr-iooor'Vo 

mi-Hr^co<HcvcTi<x)avovooo<Hrgcomc^c^Hnioior^t^vDincoor^oocwomHr^c\icococNcovomi^iHOLor'~cocM<7i!--io 

min^^'^cococococo^rTr'^^'^f'vr^Tr  ^nnirumnimnu'iininin'vr^cococoiNCMHr-1  i*  frHcjjcpinr^^rjiqiqrHcjvr^ 


r'lD'vrtNi-ioovavcococococoovO'Hcovocoi-i'vjT'OiNinr'COOH.-UNiNHHr-iocricor^vo'^rcoi-iOM^incMOvvocoovocoov 

vovovovovovoirurunimnininiDvovovovovor't^r-cococococoavoMTvooMTicriovovcocoaicocococor^r^r'r^vovovovoicuri'T 

cocococococococococorocorocococococococococococococococococococococococococococococococococococococococoroco 


888888888888888888888888888888888888888888888888888888 


r-cococococococococococococococococococococococococococococococococococococococococvavOMjvaMT'OvcvaMTicriiTicri 
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(O  O  ■'S’  VO  <-l  O  cr»  o  <£>  VO  00  o  t"  •-« CO  OJ  ■'£  t"  O 't  t"  t"  to  r>  t>  H 
t^co<NiriMc\]r-(mir)fMLOrT^'^'<rog)QinvorM"5fojLf)inrtcrit^ 
cMcor'0'TiricriOnQ^r-'X)fnvoncsoog'0>o«-(runr^covof\i 
CT>covonoonmLr)(NQQiHroLncorH^r'0'rv)voouNincOc-iTrr^ 
N(MOJC)HHOOOOOOOOQHHHH(N(NNnMn't^'f 


I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  1  I  I  I  I  I  I  I  I  I  I  I  I 


(N'£iocQa3mnQrHnr)OvDcx)r>OrHvor'^,Hmr^r^p'r)cora 

cou3ni?ioovooQOrH'<rOrHr-(vor'Oor'CJ'O^Dn<Nintir^<7iQ 

r'VDcgrHr-'>rcnor>gcQpHmLn'i)r'MOvor'(NJ2\coQ'«j,invoHLq 

cnoa'5r'0'rH(Ti'mM5'^,m<-|T3<voocgH^Dr^or>Ln>-Hinn(Nr-i 

in-^ricMfMojtNOKNrori^rinvpr^c^dtNnrsdco'p^oiH  I 

I  I  I  I  I  I  I  I  I  I  I  I  I  T  I  1  I  T  f  I  I 


ONU)rlHN 
f\l  00  r— (  *— I  rH  OJ 
r-1  fN  ID  <£>  O  !■" 

Sinoovo^o 
rncHo^t^vo 


HfNiinr'r'0'3*ooor'COCJun<NOna5in<Noo'tm 


>co<-(cp(Nfn'^,cs'r^fvojCT\onrHvo 

>CO^Qcpr^coo<N(Nvoo<Nir)r^n 


n-<rrHcavOrHcriy53'Qcpr^cbo<NCNvoo 

<N<N»X)tNfH(NnvDOinor'LnLr)iD^or^g> 

vokDvor^co<7iOiHrT«3*v£)r»^r-iroinr^c^ 


oocn  r' 
cN^inr- 


nnnnnin 


nmnrifOfo^rfrrTrTr-^r^nruriirufuDOvo'Xivo 

rifM'innnnnnnnnnnnnnnnnnn 


8888888888888888888888888888 

LnoiTiQLnoioQLnoinQirioinQinoLnQiooiOQinQinQ 

oj(nr^oojinr^OojLnr'Of^inr^O(NLnr^Of'jmr»OfNi?)r'0 

roriro^rr^Tfimmninvooo'or'r-t^r'CococoeooocT'OQ 

0\(J\O<7>(7i<7iC7l(7i<7iCri<7><7i<7>C7l(7i<7i(7iC7iC7iO<71C7iOCT\OO(7iO 
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Bearing  Error  (Deg) 

Amplitude  CdB  above  iuv/m  for  1  k W) 


Lualualei  to  San  Diego  beta-. 7. hprime-87.0 

Freq  -  23.400  kHz 

Zt  -  10.00  km  Zr  -  10.00  km 

Gamma  -  45.0  deg  Phi  *  45.0  deg 


!  HY| - !  HXi  . Bearing  Error 


Distance  (Megameters) 


Figure  4.  Plot  output  for  laterally  homogeneous  guide  sample  input. 
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arrival  error  in  degrees.  The  sample  output  shown  on  pages 
corresponds  to  the  input  which  begins  on  page  13  and  Figure 
output.  In  this  case  the  transmitter- terminator  distance  is 
Figure  2. 


30  through  37 
5  is  the  plot 
the  distance  D 


on 
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SAMPLE  OUTPUT  FOR  MOVING  TERMINATOR 


:o  oooo 

4  •  •  •  •  • 

>VOV O  VO  VO  VO 
4  CO  CO  CO  CO  CO 


C_j 

Sooooo 

Fcococococo 


orwoc^ 

Cv  r-- Cv  CO  Cvl 
fMCO^r-IH 

•  *  f  •  t 

:  i  i 


conr^rvjc 

■V'J'inr-lr 

(NrTTjTsJr 

:  *  r  *  f 


r'-vorvirsjr-' 
r-i  tn  co  rvi  •*j* 
r\jr\)coir><N 

•  •  •  •  « 
I  l  I 


O'*  T  CO  fO  «T 
COCTvr\).-(vO 
•HfNOUDfM 
•  •  •  •  • 
I  I  I 


88888 

nnnnn 


oShm't 
cr>  cv  cv  cv  cr\ 

(NtNCMCMOJ 


iiiis 

i— <  rH  i—f  rH  iH 


O  CNJC^  Qm 

t^r-cooo 

CT*  O'*  H  rH 


iXXXX 

>wwww 


_  nNcoHH 
E-l  COr-lLTlCMr-l 
_  inrHojTr'X) 

o  •  •  •  •  • 

H  '  ‘ 

co  r->  r-' r~~  vp  vo  vo 

fziM 

Q  r^nPIr- 
p  .'\i-'ininco 

+  r^CTxHCMfO 

w  •  •  •  •  • 

o  I  I 


E-t  SncouRov 

t'-cHC'J->TVO 

o  •  •  •  •  • 

H  '  ‘ 

conr'-vovpvovo 

wP??f?9 


C:  'o  co  coo 
pr-j-s-r;co 
OV>-lr-ICMn 

•  •  *  •  * 
I  I 


•P  *c 

r-poc 
coo  c 


^55533 

oroSdodd 


iWWWHH 


tT  OJrMr-IrHO 
rM  r- ii— I  rH  rH 

“  AAAAA 

Q  HfOCMt^O 

q 

O  rHrHCNCMcH 
•  «  •  «  •  • 

I  I  I  I  I 

S^J^HHgO 

8 

O  VprHHCNON 
•  vDCNjnoCn 
C\  cnc^cNr-iH 
m  •  •  •  •  • 

i  i  i  i  i 


■Vf  Cvll-HrHrHO 

.  rHr-IrHr-lrH 

“  iiiii 

O  vovocMint'- 
p  vor^t^tT'Cn 
O  (NrHCNjTTH 
•  •  •  «  •  • 
I  I  I  I  I 

8  iiiai 

O  corsir^t^co 

•  ^OOrH^1 
CV  COOVCMr-IrM 

n  •  •  •  •  • 
i  1 1 1 1 


NNHOJH 


VOHO 

\Bf\j  . 


38° 

•* 

cH  .O 

I  r*  • 

mo 

dT§ 


ooooo 


rr  o  cm  ov  vo  in  c 
ojhc 

II  II  II  II  < 


>0  O  <■»  QQ 

)•  •  o  ^  oO  * 

)ID  O  «00  **  HNCNO  ** 

«0  O  •  •  O O CO O  •  O  %» 

> fO«H  HHinOQO  •  •  •  •  m  •  -  -OO 

>  „  h  ,, 

" 11  11 11  ii  H  ii  11 11 

>Qfc  ii  e-q  ii  ii  ii  ii  ii  ii  "  "  "  ecu 


I  T ooooo! 


OvJCnJ,H<NH 


vocor^rH^o 

IDHONC^VO 

OOOCNrHCN 


< 

1  1  1  1  1 

< 

Or 

h  r->  m  cn  <n 

o 

!*fc 

CTi 

t" 

<T> 

• 

rH 

Hor'cno 

CTiHovr^o 

rH 

CVJ 

CO  CO  VO  fH  CN) 

CM 

Cm 

1  1 

U-4 

o  vdo  vo^r 
co  (NH  in  r- 
nnnncN 

«  •  •  «  • 
I  l  I  I  I 


i ncM^ovo  o 
hcjc^hi^  • 
voocNmrsj 

»  »  »  •  • 

1  «. 

3  oojco  i— i  tn  cmE 
ImnoHvf  ; 

h  o  r-  *3<  r-»  oj  t 

•  •  •  •  •  „q 
OvCTWOnO  ra 
cococococo  i— i 
w 


oojincoH 
tn^rHcon 
coco  conn 

•  •  •  #  • 
,in?  i  i  i 

3  1 

JfNncor^o 
nhnojn 
lovvocrwNr' 
•  •  •  •  * 
cvcrunr-KTv 
coco  co  cor' 
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Slab  3  R  .000  F  21.7940  A  58.500  C  39.000  M  .000  S  4.640E+00  E  81.0  T  82.0 


Ej 

£00000 

&>••••• 

OOJCMCMCMCM 

H  co  coco  co  co 


:ooooo 


)00000 
i CO  CO  CO  coco 


S  OOOOO 
•  •  •  •  • 

P CO  CO  CO  coco 

Pt^r-r-r-r- 


:ooo  00 


jvovovovovo 

<  r*  r*- r- 


S  00000 

•  •  •  •  • 


£00000 
•  •  •  •  • 

O  CM  CM  (N  CM  CM 

P  r**  r-' r- r- r- 


80rHCMCO 

corm)H 

CM  in CM  rH 


or-r-ioin 
in  r-  *3*  •sr  r— 

HnOMTUN 

r  *  f  r 


CO  COO  r-Jf-' 
^fr-HOJinO 
CO«H 


coincorHin 
c  in  <h  cmco 
r-iinOOCN 

•  •  •  •  • 

l  l  1 


in  co  in  <n  in 


Hro  r-llOLO 

in  r-  cm  «h  «h 

O  rH  lO  «H 

s  ■  r  *  r  * 


e  1 

pin^in 


ouOHnin 
in  co  co  o  cn 


SSSSS  sssss 


in  coin  cm  in 


nnHfMn 

ininnoH 

lO  CO  CO  CO  rH 


8  1 
pm^m^ 

Eoooc 


r-  cm^t  vo  co 
cm  ^Tin^rr-i 

CO  OMO  rHCM 

r #  r  r  r 


in  co^r  in  o 

OOCJriO 

niOnHn 

g  '  ‘ 

Oin^inm^" 

Kocqoo 


O  run  CTltNI 

03  in  cor-  vo 

rtOJ^NH 

1"  r  r  r  i* 


vovcminm 
vONtftrr' 
r-i  <j\  (N  t-<  r- 

g  '  '*  ’  ''  ' 

pin  coin  rv^* 
EpQOQQ 


on  cm  co  o  r- 
o  ^  in«H 
<H  CM  CO  VO  rH 


•h  iD^r 
rH  CO  C— 
O  rH  rH 
•  •  • 

K  ' 

0 10  co  in 

feQQQ 


ooogo 

imor^oQ 

r^r^r-coco 

CM  CM  CM  CM  CM 


OOOOO 

92920 

5qcd5 

CTn-lODD 

rococo  0  in 

CM  CM  CM  CM  CM 


OOOOO 

9222$ 

SecSca 

rl  CO  CM  m  CO 
Hrini-'H 
nnnnn 


OOOOO 

90900 

Soq5q 

o  vo  co  <h  p 
co  c-*  co  o  co 

CO  CO  CO  coco 


OOOOO 

90099 

c£5o5 

in  in  o  in  r- 

«—■ 1— I  <-H  O  O 

in  in  m  in  in 


o  £qpp 
in  ml-i  ovo 
vo  vo  vo  vo  in 
VO  VO  vO  vo  VO 


sss 


51-00000 

tH9<?<?Q9 


rHOUJHH 

(OiOCOOm 

cococooicn 


■'TOQOOO 

5qc  * 

(oooinm 

pamin<Ti<7i 

r~r~r'r^c^ 


jfOOOOQ 

^22222 

OClCQCl 

WOMMnrt 

PPCNVO 

iruoooo 


P'SSSQS 

illii 

(MOWHffl 
Nn'Tr'O 
rr  "d"  in 


jfOOOOO 

E_,^225j?g 

Hifliooin 
<71  C7i  pH  m  t" 
NNnnn 


inocor^oo 

OrifSl'l'P 

PJCNPaPaPa 


5POOO 

ioi 

m<Tll£> 

^><0  CO 
pH  pH  pH 


r-'oovo'xi 


VOiOvOlOlO  o  PlO^tUDU)  O  o  r-P'VOO'X)  O  o  P'P'P' 


sssss  i  sssss  *  urn  » ss 


'I'onrru) 

LOVOOCOLD 

ainn^p 

..... 
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Freq  »  21.794  kHz  Receiver  Distance  -  3821. Okm 

Zt  -  15.00  km  Zr  -  15.00  km  Gamma  »  45.0  deg  Phi  -  270.0  deg 


u 

3 

a 

c 


a> 

a 

E 

to 


0) 

u 

3 

ftp 

•n 


(M*  \  ooj  uj/Ani  aAoqe  gp)  apnqnduiv 
(63q)  J0JJ3  6u{jeag 


38 


Transmitter-Terminator  Distance  (Mm) 


IV.  RESULTS 


In  the  following  section  day  and  night  comparisons  are  given  between  calcula¬ 
tions  and  variations  in  the  DOA  at  San  Diego  measured  in  reference  2  for  the 
following  VLF/LF  transmitters:  Annapolis,  Lualualei,  Cutler,  Jim  Creek  and 
Silver  Creek. 

A.  DOA  RESULTS  FOR  ANNAPOLIS  TO  SAN  DIEGO  -  DAY 

Annapolis  transmits  at  21.4  kHz.  Geomagnetic  midpath  parameters  and  midpath 
values  for  the  ground  conductivity  and  dielectric  constant  used  for  the  cal¬ 
culations  are: 

azimuth  from  magnetic  north  =  249°, 
codip  =  23° , 

i  -*  |  -4  2 

I B |  =  magnitude  of  the  earth's  magnetic  induction  vector  -  0.537  x  10  w-/m  . 

-2 

a  =  ground  conductivity  =  10  S/M. 
e  =  earth's  dielectric  constant  =  15. 

Figure  6  shows  measured  histograms  for  two  days  in  April  1987  when  daylight 
prevailed  over  the  entire  path.  The  horizontal  scale  is  such  that  the  true 
bearing  should  be  45°  with  values  greater  than  45°  corresponding  to  signals 
arriving  from  south  of  the  true  DOA  and  values  less  than  45°  corresponding  to 
signals  arriving  from  north  of  the  true  DOA.  The  theoretical  spread,  dis¬ 
cussed  subsequently,  is  indicated  by  the  arrows  on  the  top  and  although  it 
encompasses  the  region  when  most  observations  occurred  it  does  underestimate 
the  observed  spread  indicating  that  perhaps  some  off  path  ground  or  iono¬ 
spheric  scatter  is  playing  a  role.  The  theoretical  spread  is  less  than  2°. 
To  arrive  at  the  theoretical  spread,  waveguide  calculations  were  performed  for 
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exponential  profiles  characterized  in  the  notation  of  Wait  and  Spies  (ref  10), 
by: 

t 

h  =  70  km 

/3  =  0.3  km  \  0.5  km  h  =  72  km  .  (18) 

h  =  74  km 

On  the  basis  of  a  previous  study  with  exponential  profiles  (ref.  11)  these 
profiles  appear  to  be  reasonable  bounds  for  daytime  propagation.  The  DOA  for 
these  profiles  at  the  Annapolis-San  Diego  range  (taken  to  be  3705  km)  was  then 
used  to  determine  the  spread.  Clearly  this  is  simply  one  of  many  procedures 
for  determining  the  spread  but,  unless  noted  otherwise,  it  is  the  method  used 
for  spread  comparisons  throughout  this  section. 

An  illustrative  example  of  mode  sum  behavior  as  a  function  of  range  is  shown 
in  Figure  7.  Shown  are  mode  sum  plots  for  the  scaled  transverse  and  lon¬ 
gitudinal  magnetic  intensities  H  and  H  respectively.  Also  shown  is  DOA  (or 

y  x 

.  i  < 

bearing)  error.  The  plot  is  for  the  =  0.3  km  ,  h  =  74  km  profile  and 
shows  that  even  under  daytime  conditions  large  errors  are  possible  (=  50°)  at 
ranges  <  1  Mm.  Examination  of  such  curves  for  the  six  profiles  considered 
showed  DOA  errors  less  than  about  3°  at  ranges  greater  than  about  1.5  Mm. 

B.  DOA  RESULTS  FOR  ANNAPOLIS  TO  SAN  DIEGO  -  NIGHT 

Figure  8  shows  measured  histograms  for  two  days  in  April  when  nighttime 
prevailed  over  the  entire  path.  To  arrive  at  the  theoretical  spread, 
waveguide  calculations  were  performed  for  the  following  exponential  profiles: 
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h  =  86  km 


(19) 


P  =  0.4  km'1,  0.5  km'1,  0.7  km'1-* 


87  km. 

88  km 


Again  they  were  chosen  on  the  basis  of  the  study  of  exponential  profiles  given 
in  reference  11.  The  theoretical  spread  of  about  12°  is  in  reasonable  agree¬ 
ment  with  the  observed  spreads  though  the  calculations  indicate  arrivals  from 
the  south  of  the  true  DOA  whereas  the  measurements  indicate  substantial  num¬ 
bers  of  arrivals  from  the  north. 


An  example  of  mode  sum  behavior  as  a  function  of  range  is  shown  in  Figure  9. 
The  possibility  of  DOA  errors  in  excess  of  70°  at  ranges  <  2  Mm  is  clearly 
evident.  Examination  of  mode  sum  plots  and  DOA’ s  for  the  nine  cases  con¬ 
sidered  indicates  that  errors  »  30°  are  quite  possible  at  select  ranges  in 
excess  of  about  2.5  Mm. 


C.  DOA  RESULTS  FOR  LUALUALEI  TO  SAN  DIEGO  -  DAY 

Lualualei  transmits  at  23.4  kHz.  Geomagnetic  midpath  parameters  and  midpath 
values  for  the  ground  conductivity  and  dielectric  constant  used  for  the  cal¬ 
culations  are: 

azimuth  from  magnetic  north  =  59° 
codip  -  40° 

|  B  |  =  0.413  x  10  ^w/m^ 
a  =  4  S/M 
6  =  81. 
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Figures  10  and  11  show  daytime  histograms  of  the  DOA's  measured  during  April 
of  1987.  The  true  bearing  should  be  45°  with  values  greater  than  45°  cor¬ 
responding  to  signals  arriving  from  north  of  the  true  DOA  and  values  less  than 
45°  corresponding  to  signals  arriving  from  south  of  the  true  DOA.  The 
theoretical  spread  of  about  2°  calculated  for  the  transmitter  receiver  separa¬ 
tion  of  4211  km  using  the  profiles  specified  in  Equations  (18)  is  in  quite 
good  agreement  with  the  4/10,  4/11  and  4/12  data.  The  4/1  data  indicates  a 
somewhat  larger  spread  with  the  preponderance  of  measurements  indicating 
arrivals  from  more  southerly  directions  than  predicted. 

Of  the  six  profiles  used  to  determine  the  spread  the  only  one  which  indicated 

-1  1 

an  appreciable  effect  was  the  /3  =  0.5  km  ,  h  =  74  km  case  shown  in  Figure 
12.  In  that  instance  it  will  be  seen  that  a  deep  null  in  HY  occurs  close  to 
the  receiving  site.  The  spread  indicated  on  Figures  10  and  11  could  clearly 
be  increased  by  allowing  for  a  reception  range  about  the  receiving  site  and/or 
a  more  densely  selected  set  of  profiles.  This  points  out  the  subjective 
nature  of  the  comparisons. 

D.  DOA  RESULTS  FOR  LUALUALEI  TO  SAN  DIEGO  -  NIGHT 


Figures  13  and  14  show  nighttime  histograms  of  the  DOA's  for  the  four  days  in 
April.  The  profiles  specified  in  Equations  (19)  were  used  to  determine  the 
theoretical  spreads.  The  lower  spread  was  calculated  as  described  for  the 
previous  cases.  The  small  spread  is  singular  among  all  of  the  nighttime  cases 
considered.  It  reflects  in  some  measure  the  smooth  propagation  characteris¬ 
tics  of  easterly  propagation.  If  a  loose  tolerance  of  +  200  km  in  range  about 
the  receiving  site  is  permitted  then  the  theoretical  spread  is  given  by  the 
curve  marked  "calculated  spread*".  The  reason  for  the  difference  in  the 
spread  is  evident  from  Figure  15  which  shows  the  mode  sum  and  DOA  error  for 
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the  /)  =  0 . 7  km  ,  h  =  86  km  profile.  The  error  peak  which  occurs  at  about 
4400  km  is  responsible  for  the  augmented  spread. 

E.  DOA  RESULTS  FOR  CUTLER  TO  SAN  DIEGO  -  DAY 

Cutler  transmits  at  24.0  kHz.  Geomagnetic  midpath  parameters 
values  for  the  ground  conductivity  and  dielectric  constant  used 
culations  are: 

azimuth  from  magnetic  north  =  245.5° 
codip  =  19.3° 

|  B  j  =  0.553  x  10'4  w/m2 
a  =  10'2  S/m 
€  =  15. 

Figure  16  shows  daytime  histograms  of  the  DOA's  measured  during  two  days  in 
March  1987.  The  true  bearing  should  be  45°  with  values  greater  than  45° 
corresponding  to  signals  arriving  from  south  of  the  true  DOA  and  values  less 
than  45°  corresponding  to  signals  arriving  from  north  of  the  true  DOA.  The 
theoretical  spread,  calculated  for  a  transmitter-receiver  distance  of  4454  km, 
is  somewhat  less  than  2°  and  appears  to  be  consistent  with  the  bulk  of  the 
measurements.  An  example  of  the  mode  sum  and  DOA  plots  for  this  case  is  shown 
in  Figure  17. 

F.  DOA  RESULTS  FOR  CUTLER  TO  SAN  DIEGO  -  NIGHT 

Figure  18  shows  nighttime  histograms  of  the  DOA  's  for  three  days  in  March  of 
198/.  The  profiles  specified  in  Equations  (19)  were  used  to  determine  the 
theoretical  spreads.  The  rather  large  calculated  spread  in  excess  of  30° 


and  midpath 
for  the  cal- 
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results  from  a  null  in  | HY [  close  to  the  receiving  site  (4454  km)  for  the  /?  = 
.  i  < 

0.5  km  ,  h  =  88  km  profile.  This  is  shown  in  Figure  19.  The  theoretical 
spread  is  larger  than  indicated  by  the  measurements.  The  theoretical  result 
predicts  all  arrivals  from  south  of  the  true  bearing  while  the  measurements 
show  a  substantial  number  of  arrivals  from  the  north. 

G.  DOA  RESULTS  FOR  JIM  CREEK  TO  SAN  DIEGO  -  DAY 

Jim  Creek  transmits  at  24.8  kHz.  Geomagnetic  midpath  parameters  and  midpath 
values  for  the  ground  conductivity  and  dielectric  constant  used  for  the  cal¬ 
culation  are: 

azimuth  from  magnetic  north  =  149° 
codip  =  25° 

|  B  |  =  .518  x  10  ^  w/m2 
a  =  10~2  S/M 
e  =  15. 

Figure  20  shows  daytime  histograms  of  the  DOA's  measured  during  three  days  in 
April  1987.  The  true  bearing  should  be  45°  with  values  greater  than  45° 
corresponding  to  signals  arriving  from  east  of  the  true  DOA  and  values  less 
than  45°  corresponding  to  signals  arriving  from  west  of  the  true  DOA.  The 
calculated  spread  for  a  transmitter  receiver  distance  of  1767  km  is  about  1°. 
Though  it  underestimates  the  total  measured  spread,  the  latter  is  quite  peaked 
indicating  a  spread  for  the  bulk  of  the  measurements  of  about  1°.  An  example 
of  the  mode  sum  and  DOA  plots  for  this  case  is  shown  in  Figure  21. 


H.  DOA  RESULTS  FOR  JIM  CREEK  TO  SAN  DIEGO  -  NIGHT 


Figures  22  and  23  show  nighttime  histograms  of  the  DOA's  for  five  days  in 
April  1987  and  one  day  in  December  1986.  The  profiles  specified  in  Equation 
(19)  were  used  to  determine  the  theoretical  spread  which  is  over  40°.  The 
histograms  range  from  spreads  of  about  20°  to  about  40°.  Figure  24  is  an 
example  of  the  bearing  error  behavior  with  range.  This  path  has  been  cited  as 
having  the  most  extreme  behavior  (ref.  2)  .  The  reason  for  this  is  probably 
the  relatively  short  path  length  (1767  km)  rather  than  anything  peculiar  about 
the  north  south  path.  It  is  clear  from  the  nighttime  mode  sum  plots  (e.g.  , 
Figures  9,  15,  19  and  24)  that  the  most  severe  problems,  except  for  isolated 

instances,  occur  for  ranges  <  2000  km. 

I.  DOA  RESULTS  FOR  SILVER  CREEK  TO  SAN  DIEGO  -  DAY 

Silver  Creek  transmits  at  48.5  kHz.  It  is  the  sole  LF  frequency  considered  in 
this  study.  Geomagnetic  midpath  parameters  and  midpath  values  for  the  ground 
conductivity  and  dielectric  constant  used  for  the  calculations  are: 

azimuth  from  magnetic  north  =  227.7° 
codip  =  24.9° 

|  B |  =  0.523  x  10~4w/m2 
a  •=  1.0  2s/m 
c  =  15. 


At  this  LF  frequency  the  study  reported  in  reference  11  indicates  that  the 
following  daytime  profiles  cover  the  range  of  interest 
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Figures  25,  26  and  27  show  daytime  histograms  for  six  days  of  measurement  in 
November  1986  and  one  day  in  December  1986.  Values  greater  than  45°  cor¬ 
respond  to  signals  arriving  from  south  of  the  true  bearing  and  values  less 
than  45°  correspond  to  signals  arriving  from  north  of  the  true  bearing.  The 
theoretical  spread  of  about  2°  calculated  on  the  basis  of  the  above  profiles 
at  a  range  of  1989  km  is  consistent  with  the  observations.  The  calculated 
result  indicates  arrivals  from  the  north  whereas  the  measurements  indicate 
arrivals  from  both  north  and  south  of  the  true  bearing.  A  sample  mode  sum  and 
DOA  plot  for  this  case  is  given  in  Figure  28. 


J.  DOA  RESULTS  FOR  SILVER  CREEK  TO  SAN  DIEGO  -  NIGHT 


Figures  29,  30  and  31  show  nighttime  histograms  for  the  seven  days  of  measure¬ 
ments.  Under  nighttime  conditions  at  this  LF  frequency  the  following  six 
profiles  have  been  used. 

h  =  88  km. (21) 

Again  thev  were  chosen  on  the  basis  of  the  study  of  exponential  profiles  given 
in  reference  11. 

The  theoretical  spread  based  on  these  profiles  indicates  a  preponderance  of 
arrivals  from  the  north  with  the  spread  being  over  50°.  The  measurements 
indicate  spreads  ranging  between  about  20°  and  30°  with  substantial  arrivals 
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from  both  the  North  and  South.  The  reason  for  the  large  calculated  spread  is 
the  deep  null  in  | HY |  at  about  2  Mm  shown  in  Figure  32.  Failure  of  the 
measurements  to  show  such  an  effect  could  indicate  that  either  the  profiles 
were  inadequate  or  that  too  few  modes  were  used.  With  respect  to  the  latter, 
in  all  cases  the  modal  search  rectangle  was  60°  to  90°  for  the  real  part  of 
the  eigenangle  and  0  to  -4°  for  the  imaginary  part.  Examination  of  the  mode 
set  in  this  instance  indicates  that  the  dominant  mode  would  be  at  least  30  dB 
higher  than  the  first  neglected  mode  at  the  receiving  site  so  that  profile 
inadequacy  seems  the  more  likely  explanation. 
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V.  CONCLUSIONS 


A  modification  of  an  earlier  program  has  been  used  to  calculate  DOa  errors  due 
to  geomagnetic  field  effects  at  San  Diego  for  the  following  VLF/LF  transmit¬ 
ters:  Annapolis  (21.4  kHz),  Lualualei  (23.4  kHz),  Cutler  (24  kHz),  Jim  Creek 

(24.8  kHz)  and  Silver  Creek  (48.5  kHz).  Calculated  spreads  of  the  DOA  have 
been  based  on  a  selection  of  best  fit  exponential  profiles  determined  in  a 
previous  study  (ref.  11).  This  is  subjective  in  the  sense  that  it  is  just  one 
of  many  procedures  that  could  be  used  to  determine  the  DOA  spread.  The  cal¬ 
culations  have  been  compared  with  corresponding  variations  in  the  DOA  reported 
in  reference  2.  Calculated  daytime  spreads  were  about  2°.  Measured  daytime 
spreads  were  somewhat  larger  than  this,  though  the  bulk  of  the  measurements 
usually  is  confined  to  a  range  of  several  degrees.  In  agreement  with  measure¬ 
ments  and  expectations,  the  nighttime  calculated  DOA  spreads  were  much  greater 
than  the  daytime  spreads.  A  quantitative  comparison  between  calculated  and 
measured  nighttime  spreads  is  shown  below: 


Calculated 
Nighttime  Spread 


Measured 

Nighttime  Spread 


Annapolis  =  11° 
L-ialualei  ~  2° 
Cutler  =:  32° 
Jim  Creek  ~  42° 
Silver  Creek  =  53° 


=  25° 


~  20° 


=  22°  to  ~  42° 
-  33°  to  47° 


In  addition,  the  nighttime  calculations  for  Annapolis  and  Cutler  yield  an 
excessive  number  (relative  to  measurement)  of  arrivals  from  the  south,  the 
calculations  for  Jim  Creek  indicate  an  excessive  number  from  the  east  and  the 
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calculations  for  Silver  Creek  indicate  an  excessive  number  of  arrivals  from 
the  north. 

Generally  there  is  no  obvious  trend  in  DOA  variation  with  frequency  in  the  VLF 
band,  though  on  the  basis  of  the  calculations  the  nighttime  DOA  variation  is 
expected  to  be  more  severe  for  the  LF  frequency  than  for  the  VLf  frequencies 
at  ranges  >  5000  km.  Calculations  and  measurements  indicate  that  the 
nighttime  westerly  paths  are  subject  to  more  DOA  variation  than  the  single 
easterly  path  considered.  This  is  believed  to  be  due  to  geomagnetic  field 
effects  and  not  to  conductivity  differences  for  the  easterly  and  westerly 
paths  (ref.  12).  Reference  2  singles  out  Jim  Creek  as  a  path  of  excessive  DOA 
variation.  On  the  basis  of  the  calculations  it  is  believed  that  this  is 
because  of  the  relatively  short  path  (1767  km)  and  not  because  of 
peculiarities  associated  with  north-south  propagation. 
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Figure  6.  Histograms  for  Annapolis  to  San  Diego  -  Day. 
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Figure  8.  Histograms  for  Annapolis  to  San  Diego  -  Night. 
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Figure  9.  Sample  range  plot  for  Annapolis  to  San  Diego  -  Night. 
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Figure  10.  Histograms  for  Lualualei  to  San  Diego  -  Day. 
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Figure  12.  Sample  range  plot  for  Lualualei  to  San  Diego  -  Day. 
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Figure  13.  Histograms  for  Lualualei  to  San  Diego  -  Night, 
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Figure  14.  Histograms  for  Lualualei  to  San  Diego  -  Night. 


60 


Bearing  Error  (Deg) 

Amplitude  (d9  above  iuv/m  for  1  kW) 


Lualualei  to  San  Diego  beta-. 7,  hprime=B6.0 

Freq  -  23.400  kHz 

Zt  -  0.00  km  Z r  -  0.00  km 

Gamma  -  0.0  deg  Phi  -  0.0  deg 


0  2  4  6  8  10 

Distance  (Megameters) 


Figure  15.  Sample  range  plot  for  Lualualei  to  San  Diego  -  Night. 
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Figure  17.  Sample  range  plot  for  Cutler  to  San  Diego  -  Day. 
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Figure  19.  Sample  range  plot  for  Cutler  to  San  Diego  -  Night. 
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Figure  24.  Sample  range  plot  for  Jim  Creek  to  San  Diego  -  Night. 
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Figure  25.  Histograms  for  Silver  Creek  to  San  Diego  -  Day. 
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Figure  26.  Histograms  for  Silver  Creek  to  San  Diego  -  Day. 
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Figure  27.  Histograms  for  Silver  Creek  to  San  Diego  -  Day. 
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Figure  28.  Sample  range  plot  for  Silver  Creek  to  San  Diego  -  Day. 
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Histograms  for  Silver  Creek  to  San  Diego  -  Night 
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Figure  30.  Histograms  for  Silver  Creek  to  San  Diego  -  Night. 
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Figure  31.  Histograms  for  Silver  Creek  to  San  Diego  -  Night 
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Figure  32.  Sample  range  plot  for  Silver  Creek  to  San  Diego  -  Night. 
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APPENDIX  A:  PROGRAM  LISTING 


parameter  (maxmds— 30) 

IMPLICIT  REAL  *8(A-H,0-Z) 

COMMON/HGTEMP/FF1 (25, maxmds ) , FF3 (25 , maxmds ) 

COMMON/TERM/NT , NTR 

COMMON/CAP/CAP 1(25, maxmds , maxmds ) , TNORM (25, maxmds , maxmds ) 
COMMON/MCINPT/THETA(25 , maxmds) , FOFR(25 , maxmds) , 

XTRA (3,3,25, maxmds ) ,TOPHT(25) , 

XVAL(25) , FREQ , RHOMAX , RHOMIN , DELRHO , DELTAX, epsr (25) , 
S IGMA ( 2 5 ) , NRSLAB , NRMODE ( 2 5 ) , NTMAX 
COMMON/MCSTOR/A(25 .maxmds .maxmds) , S(25 .maxmds) ,C(25 .maxmds) , 

$  NTHSQ (25), KVRAOT , KVRATT , 

$  AVRKOT , AVRKTT , CONST , OMEGA , WAVENO 

COMMON /MC  PLOT /R (400) ,DBy(400) ,DBx(400) ,ANGchi(400) , 

$  IDPLOT.ISUB 

COMMON/XPLOT /XMIN , xmax ,Xtic , YMIN , ymax , Ytic , SIZEX , SIZEY 

COMMON /HG I NPT /GAMMA , PHI , ZT , ZR , SINGAM, COSGAM , SINPHI , COSPHI 

common/pltflg/iplf lg 

DIMENSION  Z(2) 

character*^  bcd(20) 

character*68  idplot 

REAL*4  R, DBy , dbx , angchi 

REAL* 4  XMIN, xmax, Xtic, YMIN, ymax, Ytic, SIZEX, SIZEY 
COMPLEX* 16  FF1.FF3 

COMPLEX*16  THETA , A , S , C , FOFR , IM , CAPI , TNORM 
complex*16  temp 

COMPLEX*16  XTRA, NTHSQ, T1,T2,T3,T4 
REAL*8  KVRAOT , KVRATT 
data  IM/(0 . 0D0 , 1 . 0D0) / 

NAMELIST/DATUM/  RHOMAX, 

&  RHOMIN, DELTAX, NRSLAB, NTMAX ,XVAL, 

$  DELRHO, IFIRST, LAST, IPLTOP , iplf lg , 

$  XMIN, xmax, Xtic, YMIN, ymax, Ytic .SIZEX, SIZEY, 

$  GAMMA , PHI , ZT , ZR , IPRNTA 
$  , INTFLG 

DATA  TWOPI/6 . 283185D0/ , VELITE/2 . 997928D5/ , ALPHA/3 . 14D-4/, 

$  DEGRAD/1. 745329D-2/ 

DATA  LAST/0/, I PRNTA/0/ 

C 

INTFLG=0 
iplflg-0 
10  continue 

READ (5, DATUM) 
print  datum 
write(*. '(/)') 

GAMMA  -  GAMMA*DEGRAD 
PHI  -  PHI*DEGRAD 
SINGAM  -  DS IN (GAMMA) 

COSGAM  -  DCOS (GAMMA) 

SINPHI  -  DS IN ( PHI ) 

COSPHI  -  DCOS (PHI) 
read( 5 , ' (a68)  ' )  idplot 
write(*, ' (lx,a68 ,/) ' )  idplot 
DO  M-l, NRSLAB 
nrmode (m)=l 

re ad (5, 1020)  rr, f  f , aa , cc , bb , ss.ee, topht ' m) 
print  1022,m,rr,ff,aa,cc,bb,ss,ee, topht (m) 

PRINT  102 

if(m  . gt .  1  .and.  ff  .ne.  freq  )  then 
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print  *, 'Check  input  records  -  mis-match  on  freq’ 
stop 
endif 
freq-ff 
sigma(m)-ss 
epsr (m)-ee 
205  k-nrmode (m) 

read(5 , 1023)  indxl , temp , ipoll , tl , t2 
if(indxl  .gt.  0)  then 

read(5 , 1023)  indx2 , theta(m,k) , ipol2 , t3 , t4 
c  get  ey/hy 

if(ipoll  .eq.  1)  then 
fofr(m,k)-t3/tl 
else 

fofr(m,k)=t2/( t3*t4) 
end  if 

PRINT  103 ,THETA(M,K) , Tl , T2 , T3 , T4 , F0FR(M , K) ,T0PHT(M) 
S(M,K)  -  zsin(THETA(M,K)*DEGRAD) 

C(M,K)  -  zcos (THETA (M,K)*DEGRAD) 

XTRA(1 , 1 ,M,K)  -  T1*S (M , K)**2 
XTRA(1 , 2 ,M,K)  -  T1*S(M,K) 

XTRA(1,3,M,K)  -  -T3*S(M,K) 

XTRA(2,1,M,K)  -  -T1*S(M,K) 

XTRA(2 , 2 ,M, K)  -  -Tl 

XTRA(2 , 3 ,M,K)  -  T3 

XTRA(3 , 1 ,M,K)  -  -T3*T4*S(M,K) 

XTRA( 3 , 2 , M , K)  -  -T3*T4 
XTRA(3 , 3 ,M,K)  -  T2 
nrmode (m) -nrmode (m) +1 
go  to  205 
endif 

nrmode (m) -nrmode (m) - 1 
enddo 
NT  -  1 

WAVE  NO  -  TWO  PI*1000 . 0*FREQ/VELITE 
CONST  -  0 . 03248*WAVE  NO/DSQRT(FREQ) 

OMEGA  -  TWO  PI*FREQ*1000. 

KVRAOT  -  DEXP ( DLOG (WAVE  NO/ALPHA) /3 . ) 

KVRATT  -  KVRAOT** 2 
AVRKOT  -  1. /KVRAOT 
AVRKTT  -  AVRKOT** 2 *0 . 5 
DO  130  L=>1 ,  NRSLAB 
130  NTHSQ(L)  -  1 . +ALPHA*TOPHT(L) 

IFLG  -  0 

DO  135  M-l, NRSLAB 
IF(M  .NE.  1)  IFLG-1 

135  CALL  HTINTL(CAPI ,TNORM, IFLG ,M, INTFLG) 

Z(l)  -  ZT 
Z ( 2 )  -  ZR 
CALL  HTGAIN(Z) 

DO  136  M-l, NRSLAB 
DO  136  K-l, NRMODE (m) 

XTRA( 1 , 1 , M , K)  -  XTRA( 1 ,1,M,K)*FF1(M,K)**2 

XTRA( 1 , 2 , M , K)  -  XTRA( 1 , 2 ,M , K)*FF1 (M, K)**2 

XTRA(1,3,M,K)  -  XTRA(1,3,M,K)*FF1(M,K)*FF3(M,K)/F0FR(M,K) 

XTRA( 2 , 1 , M , K)  -  XTRA(2 , 1 , M , K) *FF1 (M , K) **2 

XTRA(2 ,2 ,M,K)  -  XTRA( 2 , 2 , M , K) *FF1 (M , K) **2 

XTRi\(  2 , 3  ,  M  ,  K)  -  XTRA(  2 , 3  ,  M ,  K)  *Fc'l  (M,K)*FF3(M,K)/FOFR(M,K) 
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XTRA(3,1,M,K)  -  XTRA(3,1,M,K)*FF1(M,K)*FF3(M,K)/F0FR(M,K) 
XTRA(3 , 2 ,M,K)  -  XTRA(3 . 2 ,M,K)*FF1(M,K)*FF3(M,K)/F0FR(M,K) 
XTRA( 3 , 3 ,M,K)  -  XTRA(3 , 3 ,M,K)*FF3 (M,K)**2/(FOFR(M,K)**2) 
136  CONTINUE 

if(nrslab  .eq.  1)  then 
do  j-1 .nrmode (1) 
do  k-1 , nrmode (1) 
a(l,k,j)-(0. ,0) 

if(k  .eq.  j)  a(l,k,j)-(l. 0,0.0) 
enddo 
enddo 
ntr-1 
go  to  405 
endif 

118  IF(XVAL(2)  .GE.  0.)  GO  TO  111 
DO  112  L-l, NRSLAB 
DO  112  J-l.NRMODE(l) 

DO  112  K-l.NRMODE(l) 

112  A(L,K, J)  -  0.0 

DO  113  L— 2 .NRSLAB 

IF (XVAL(L)  .GE.  0.)  GO  TO  114 

113  CONTINUE 
NTR  -  nrslab 
GO  TO  117 

114  NTR  -  L-l 
117  CONTINUE 

DO  116  J-1, NRMODE (ntr) 

DO  116  K-l.NRMODE(ntr) 

116  IF(K  .EQ.  J)  A(NTR,K,J)-(1. 0,0.0) 

DO  401  M-ntr, nrslab 
401  CALL  MCSTEP(M) 

IF(IPRNTA  .EQ.  0)  GO  TO  91 
PRINT  905 

DO  451  L-ntr, nrslab 

PRINT  900 ,L 

DO  451  J-1, NRMODE (1) 

DO  451  K-1, NRMODE (1) 

PRINT  901 , J ,K,A(L, J ,K) 

451  CONTINUE 

91  IF(IPLTOP  .EQ.  1)  CALL  MCFLD 
IF ( IPLTOP  .EQ.  2)  CALL  MCFLD 2 
NT  -  NT+1 

DO  106  ME=2, nrslab 
106  XVAL(ME)  -  XVAL(ME)+DELTAX 

IF(XVAL(NTR)  .GE.  0.  .AND.  NT  .LE.  NTMAX)  GO  TO  118 
IF(NT  .LE.  NTMAX)  GO  TO  91 
IF(LAST  .EQ.  0)  GO  TO  10 
stop 

111  NTR  -  1 

DO  120  L  -  1, NRSLAB 
DO  120  K  -  1,  NRMODE(l) 

DO  120  J  -  1,  NRMODE(l) 

A(L,K, J)  -  0.0 

IF( 1  .eq.  1  .and.  K  .EQ.  J)  A(l,j,k)  -  (1. 0,0.0) 

120  CONTINUE 

C  THE  LOOP  400  DETERMINES (A) 

C  IN  SUCCESSIVE  SLABS. 

DO  400  M  -  1, NRSLAB 
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CALL  MCSTEP(M) 

IF(NRSLAB  .LE.  1)  stop 
400  CONTINUE 

405  continue 

IF(IPRNTA  .EQ.  0)  GO  TO  90 

PRINT  905 

DO  450  L-l , NRSLAB 

PRINT  900, L 

DO  450  J-l.NRMODE(l) 

DO  450  K— 1 , NRMODE(l) 

PRINT  901 , J , K, A(L, J , K) 

450  CONTINUE 

90  IF( IPLTOP  .EQ.  1)  CALL  MCFLD 
IF( IPLTOP  .EQ.  2)  CALL  MCFLD 2 
NT  -  NT  +  1 
DO  105  ME  -  2 , nrslab 
XVAL(ME)  -  XVAL(ME)  +  DELTAX 
105  CONTINUE 

IF(NT  .LE.  NTMAX)  GO  TO  90 
IF(LAST  .EQ.  0)  GO  TO  10 
stop 

102  FORMAT ( 5X , 'THETA' ,15X, 'Tl' ,20X, 'T2' ,20X, 'T3' ,20X, 'T4' ,20X, ' FOFR' , 

$  lOX.'TOPHT') 

103  FORMAT ( '  '.2F7.3.2X,  5 (2D10 . 3 , 2X) , F4 . 1) 

201  FORMAT (20A4) 

202  FORMAT ('  '.20A4) 

900  FORMAT (1H  ,14X, 

$  'A  -  TOTAL  CONVERSION  COEFFICIENTS ', 6X ,' SLAB  NUMBER  -  ',12,/) 

901  FORMAT (14X, '  J-',I2,5X,'  K  -' ,  12, 5X,'  A-' , (E15.5.E15.5) ,/) 

905  FORMAT (1H1) 

1020  format (lx ,f7.0,3(2x,f8.0) , 2 (2x , elO . 0) , 2 (2x , e5 . 0) ) 

1022  format(/,'  Slab  ',12,'  R',f7.3,'  F',f8.4,'  A’,f8.3,'  C',f3.3,’  M' , 
$  f 6 . 3 , '  S '  ,  IpelO . 3 , '  E ' , Opf 5 . 1 , '  T',f5.1) 

1023  format( il , 2f9 . 0 , il , 4el5 . 0) 

END 
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SUBROUTINE  HTINTL(CAPI .NORM, IFLG.M, INTFLG) 

C  CALCULATE  NORMALIZATION  INTEGRALS  AND  INTEGRALS  OF  HEIGHT  GAINS  IN 
C  ADJACENT  SLABS . 

parameter  (maxmds-30) 

IMPLICIT  REAL  *8(A-H,0-Z) 

C0MM0N/MCINPT/THETA(25 .maxmds) , FOFR(25 .maxmds) , 

$  XTRA(3 ,3,25 .maxmds) ,TOPHT(25) , 

$  XVAL(25) , FREQ , RHOMAX , RHOMIN , DELRHO , DELTAX , epsr (25 ) , 

$  SIGMA(25) , NRSLAB ,NRMODE(25) ,NTMAX 

COMMON/MCSTOR/A(25 , maxmds .maxmds) , S (25 .maxmds) , C ( 2  5 , maxmds) , 

$  NTHSQ(25) , KVRAOT , KVRATT , 

$  AVRKOT . AVRKTT , CONST, OMEGA. WAVENO 

COMPLEX* 16  NTHSQ 

COMPLEX*16  PTHA ,H1TA,H2TA,H1PRTA,H2 PRTA , HYTHA ( maxmds ) , 

$  EYTHA(maxmds) , HYTHPA(maxmds ) , EYTHPA( maxmds) 

COMPLEX*16  THETA , FOFR , A, S ,C , SSQ, CSQ, IM, NGSQ , 

&  SQROOT , RTIORT , PO , PTH , H10 , H20 , H1PRM0 , H2PRM0 , CAPH10 , CAPH20 , 

&  AlST , A2ND , A3RD , A4TH , DEN12 , DEN34 , DENMF , NURMF , 

&  HIT, H2T.H1PRMT.H2PRMT.HYTH (maxmds) , EYTH(maxmds) .HYTHPR (maxmds) , 
$  EYTHPR( maxmds) , 

&  HYOPR(maxmds) , EYOPR(maxmds ) , EYO (maxmds ) ,MULT, FAC1 , FAC2 , 

$  NORM(25 .maxmds .maxmds) , PS (maxmds) , 

$  CAPI (25 .maxmds .maxmds) , PHYTH( maxmds ) , PHYTHP( maxmds) , 

$  PEYTH (maxmds) , PEYTHP (maxmds) , PEYO (maxmds) , 

&  PEYO  PR ( maxmds ) , PHYOPR ( maxmds ) , XTRA 

REAL* 8  KVRAOT , KVRATT 
DATA  EPSLNO/8 .85434D-12/ 
data  IM/(0 . DO . 1 . DO) / 

C 

C 

DO  100  K  -  1.  NRMODE(m) 

SSQ  -  S(M,K)**2 
CSQ  -  C (M , K)**2 

c  NGSQ  -  (EPSLON(M)  -  IM*SIGMA(M)/OMEGA)/EPSLNO 

ngsq-epsr(m) - (im*sigma(m) )/(omega*epslnO) 

SQROOT  -  zsqrt (NGSQ  -  SSQ) 

RSQR  -  SQROOT 

IF(RSQR  .LT.  0.)  SQROOT-- SQROOT 

RTIORT  -  1 . /NGSQ*SQROOT 

PO  -  KVRATT*CSQ 

PTH  -  KVRATT* ( NTHSQ (M) - SSQ) 

CALL  MDHNKL(P0 . H10 , H20 , HlPRMO , H2PRM0) 

CAPH10  -  HlPRMO  +  AVRKTT*H10 

CAPH20  -  H2PRM0  +  AVRKTT*H20 

AlST  -  CAPH20  -  IM*RTIORT*KVRAOT*H20 

A2ND  -  CAPH10  -  IM*RTIORT*KVRAOT*H10 

A3RD  -  H2PRM0  -  IM*KVRAOT*SQROOT*H20 

A4TH  -  HlPRMO  -  IM*KVRAOT*SQROOT*H10 

DEN12  =  H20*A2ND  -  H10*A1ST 

DEN34  -  H20*A4TH  -  H10*A3RD 

CALL  MDHNKL( PTH , HIT . H2T , H1PRMT , H2PRMT) 

HYTH(K)  -  (H2T*A2ND  -  H1T*A1ST)/DEN12 
EYTH(K)  -  (H2T*A4TH  -  H1T*A3RD)/DEN34*F0FR(M , K) 

HYTHPR(K)  =  (H2PRMT*A2ND  -  H1PRMT*A1ST)/DEN12 
EYTHPR(K)  -  (H2PRMT*A4TH  -  H1PRMT*A3RD)/DEN34*F0FR(M , K) 

HYOPR(K)  -  (H2PRM0*A2ND  -  H1PRM0*A1ST)/DEN12 
EYOPR(K)  -  (H2PRM0*A4TH  -  H1PRM0*A3RD)/DEN34*F0FR(M,K) 

I F ( IFLG  .EQ.  0)  GO  TO  100 
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PTHA  -  KVRATT*(NTHSQ(M-1) -SSQ) 

CALL  MDHNKL(PTHA , H1TA , H2TA , H1PRTA , H2PRTA) 

HYTHA(K)  -  (H2TA*A2ND-H1TA*A1ST) /DEN12 
EYTHA(K)  -  (H2TA*A4TH-H1TA*A3RD)/DEN34*F0FR(M,K) 

HYTHPA(K)  -  (H2PRTA*A2ND-H1PRTA  *A1ST)/DEN12 
EYTHPA(K)  -  (H2PRTA*A4TH-H1PRTA*A3RD)/DEN34*F0FR(M,K) 

100  EYO(K)  -  FOFR(M,K) 

IF(INTFLG  .EQ.  1)  PRINT  906,  M 
DO  240  J  -  1 ,NRMODE(m) 

DO  240  K  -  1 , NRMODE(m) 

IF( J  .EQ.  K)  GO  TO  120 

MULT  -  AVRKOT / ( ( S (M , J )  -  S (M,K) )*WAVENO) 

FAC1  -  EYTH(K)*EYTHPR(J)  -  EYTH ( J ) *EYTHPR (K)  +  HYTH ( K ) *HYTHPR ( J ) 

$  -HYTH(J)*HYTHPR(K) 

FAC 2  -  -EYO (K)*EY0PR(J)  +  EY0(J)*EY0PR(K)  -  HYOPR(J)  +  HYOPR(K) 
NORM(M, J ,K)  -  MULT* ( FAC 1+FAC2) 

IF(INTFLG  .EQ.  1)  PRINT  908 ,M, J ,K,NORM(M, J ,K) 

GO  TO  240 

120  MULT  -  2.0*S(M,J)*KVRA0T/WAVEN0 
PTH  -  KVRATT*(NTHSQ(M) -S(M, J)**2) 

PO  -  KVRATT*C (M , J ) **2 

FAC1  -  EYTHPR(J)**2  +  HYTHPR(J)**2  +  PTH*(EYTH(J)**2  +  HYTH(J)**2) 
FAC 2  -  -EY0PR(J)**2  -  HY0PR(J)**2  -  P0*(EY0(J)**2  +  1.0) 
NORM(M,J,K)  -  MULT* ( FAC 1+FAC2) 

IF(INTFLG  .EQ.  1)  PRINT  908 ,M, J ,K,NORM(M, J ,K) 

240  CONTINUE 

IF  (IFLG  .EQ.  0)  GO  TO  500 
DO  400  K  -  1,  NRMODE(ra) 

DO  400  J  -  1,  NRMODE(m-l) 

MULT  -  AVRKOT/ ((PS(J)  -  S(M,K) )*WAVENO) 

FAC1  -  EYTHA ( K ) *PEYTHP ( J ) -  PE YTH (J) *EYTH PA ( K ) 

$+HYTHA (K) *PHYTHP ( J ) - PHYTH ( J ) *HYTHPA(K) 

FAC2  -  -EYO (K)*PEY0PR(J)  +  PEY0(J)*EY0PR(K)  -PHYOPR(J)  +  HYOPR(K) 
CAPI (M,K, J)  -  MULT*(FAC1+FAC2) 

IF(INTFLG  .EQ.  1)  PRINT  910 ,M ,K ,J , CAPI (M, K, J) 

400  CONTINUE 

500  DO  600  J  -  1,  NRMODE(ra) 

PS(J)  -  S (M , J ) 

PHYTH(J)  -  HYTH(J) 

PHYTHP(J)  =  HYTHPR(J) 

PEYTH(J)  -  EYTH(J) 

PEYTHP(J)  =  EYTHPR(J) 

PHYOPR(J)  -  HYOPR(J) 

PEiO(J)  -  EYO(J) 

PEYOPR(J)  =  EYOPR(J) 

600  CONTINUE 
720  CONTINUE 
RETURN 

906  FORM*.T(’0' ,20X, 'INTEGRALS  IN  SLAB' ,  13,/) 

908  F0RMAT(21X, 'NORM(' ,11, ' , '  ,11, ' , ' ,11, ')  -'.2D13.6) 

910  FORMAT (2 IX, ' CAPI (',11, ',',11, ',',11,')  -'.2D13.6) 

END 
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SUBROUTINE  HTGAIN(Z) 

C  COMPUTE  EZ , EX , EY  HEIGHT  GAINS  FOR  TRANSMITTER  AND  RECEIVER, 
parameter  (maxmds-30) 

IMPLICIT  COMPLEX*16 (A-H.O-Z) 

COMMON/HGTEMP /FF1 (25 , maxmds ) , FF3 (25 , maxmds ) 

COMMON /HTGN /F ( 4 , 25 , maxmds , 2 ) 

COMMON/MCINPT/THETA(25,,naxmds)  ,FOFR(25, maxmds) , 

XTRA(3 ,3,25, maxmds) ,TOPHT(25) , 

XVAL(25) , FREQ , RHOMAX , RHOMIN , DELRHO , DELTAX ,epsr(25) , 
SIGMA(25) ,NRSLAB,NRMODE(25) ,NTMAX 
COMMON/MCSTOR/A(25 .maxmds .maxmds) , S (25 .maxmds) , C(25 .maxmds) , 

$  NTHSQ (25), KVRAOT , KVRATT , 

$  AVRKOT.AVRKTT, CONST, OMEGA, WAVENO 

COMPLEX* 16  zsqrt 
REAL* 8  DEXP 

COMPLEX*16  NGSQ.IM, NTHSQ 

REAL* 8  XVAL , FREQ , RHOMAX , RHOMIN , DELRHO , DELTAX , epsr .SIGMA 

REAL* 8  KVRAOT , KVRATT , AVRKOT , AVRKTT , CONST , OMEGA , WAVENO 

REAL* 8  Z(2) .EPSLNO, ALPHA, FAC1 

REAL* 8  RSQR 

REAL* 8  TOPHT 

data  IM/ (0 . DO , 1 . DO) / 

data  EPSLNO/8 . 85434D- 12/ , ALPHA/3 . 14D-4/ 
c 

DO  100  M-l.NRSLAB 

c  NGSQ  -  (EPSLON (M) - IM*SIGMA(M) /OMEGA) /EPSLNO 
ngsq-epsr(m) - (im*sigma(m) ) /(oraega*epslnO) 

DO  100  K-l.NRMODE(m) 

SSQ  -  S (M,K)**2 
SQROOT  -  zsqrt (NGSQ- SSQ) 

CSQ  -  C(M,K)**2 
RSQR  -  SQROOT 

IF (RSQR  .LT.  0.)  SQROOT=- SQROOT 
DO  100  IZ-1,2 

Q  -  KVRATT* (CSQ+ALPHA*Z(IZ)) 

QO  -  KVRATT*CSQ 

CALL  MDHNKL(Q0 , H10 , H20 , H1PRM0 , H2PRM0) 

CALL  MDHNKL(Q , HI , H2 , H1PRM , H2PRM) 

CAPH10  -  H1PRM0+AVRKTT*H10 
CAPH20  -  H2PRM0+AVRKTT*H20 
FAC 2  -  IM*KVRAOT*SQROOT 
FAC3  -  FAC2/NGSQ 
FI  -  - (CAPH20-FAC3*H20) 

F2  -  CAPH10 - FAC3*H10 
F3  -  - (H2PRM0-FAC2*H20) 

F4  -  H1PRM0 - FAC2*H10 

FAC1  -  DEXP(ALPHA/2 . *Z ( IZ) ) 

F(1,M,K,IZ)  -  FAC1* ( F1*H1+F2*H2 ) 

F(2,M,K, IZ)  -  ALPHA/ ( IM*2 . *WAVENO)*F (1,M,K,IZ)+1. /IM*AVRK0T*FAC1* 
$  ( F1*H1PRM+F2*H2PRM) 

F(3,M,K,IZ)  =  F3*H1+F4*H2 

f(4,m,k,iz)  =  - im*avrkot*(f3*hlprm+f4*h2prm) 

FFl(M.K)  =  F1*H10+F2*H20 
FF3 (M , K)  -  F3*H10+F4*H20 
F( 1 ,M,K , IZ)  -  F(1,M,K,IZ)/FF1(M,K) 

F(2 ,M,K. IZ)  -  F( 2 , M , K , IZ) /FF1 (M , K) 

F( 3 , M , K , IZ)  -  F(3,M,K,IZ)*FOFR(M,K)/FF3(M,K) 

F(4.M,K,IZ)  -  F(4,M,K,IZ)*FOFR(M,K)/FF3(M,K) 
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100  CONTINUE 
RETURN 
END 
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SUBROUTINE  MCSTEP(M) 

C  COMPUTE  GENERALIZED  MODE  CONVERSION  COEFFICIENTS, 
parameter  (maxmds-30) 

IMPLICIT  REAL  *8(A-H,0-Z) 

COMMON/TERM/NT , NTR 

COMMON/CAP/CAPI (25, maxmds , maxmds ) , TNORM (25, maxmds , maxmds ) 
COMMON/MCINPT/THETA(25 .maxmds) , FOFR( 25 .maxmds) , 

XTRA(3 ,3,25 .maxmds) ,TOPHT(25) , 

XVAL(25) , FREQ ,  RHOMAX , RHOMIN , DELRHO , DELTAX ,epsr(25) , 
SIGMA(25) , NRSLAB,NRMODE(25) , NTMAX 
COMMON/MCSTOR/A(25 , maxmds .maxmds) , S (25 .maxmds) , C(25 .maxmds) , 

$  NTHSQ(25) , KVRAOT , KVRATT , 

$  AVRKOT.AVRKTT, CONST, OMEGA .WAVENO 

COMPLEX*16  zexp 
COMPLEX* 16  NTH SQ 

COMPLEX*16  THETA , FOFR , A , S , C , TNORM , CAPI , 

$  IM, B (maxmds) ,ANS (maxmds) ,TS (maxmds .maxmds)  XTRA 
REAL*8  KVRAOT , KVRATT 
REAL*4  ERR 
data  IM/ ( 0 . DO , 1 . DO ) / 

C 

C 

MP  -  M-l 

IF(M  -EQ.  NTR)  RETURN 
DO  N-l, maxmds 
B(N)  -  (0 . , 0 . ) 
enddo 

IF (MP  .ne.  NTR)  then 
DO  K-l  .NRMODE(ntr) 

DO  L-l .NRMODE(m) 

DO  J-l, NRMODE (m) 

TS(L,J)  -  TNORM(M.L.J) 
enddo 
enddo 

do  1-1 .nrmode(m) 
do  j— 1 , nrmode (mp) 

B(L)— B (L)+A(MP , J . K)*zexp( - IM*WAVEN0*S (MP , J )* (XVAL(M) - 
&  XVAL(MP) ) )*CAPI (M , L, J) 

enddo 
enddo 

CALL  CLINEQ(TS,B,ANS, NRMODE (m) .maxmds , 0 , ERR) 

DO  I-l.NRMODE(m) 

c  A(M, I ,K)  -  ANS (I )*S (M , I ) /S (NTR , K) 

A(M.I.K)  -  ANS ( I ) 
enddo 

DO  N-l, maxmds 
B(N)  -  (0 . , 0 . ) 
enddo 
enddo 
else 

DO  K-l , NRMODE(mp) 

DO  LL-1, NRMODE (m) 

DO  L-l.NRMODE(m) 

TS(LL,L)  -  TNORM (M.LL.L) 

B(L)  -  CAPI(M.L.K) 
enddo 
enddo 

CALL  CLINEQ(TS , B , ANS , NRMODE (m) , maxmds , 0 , ERR) 
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DO  J-l,NRMODE(m) 

A(M, J,K)  -  ANS(J)*S(M, J)/S(NTR,K) 
A(M,J,K)  -  ANS(J) 
enddo 
enddo 
endif 
RETURN 
END 
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SUBROUTINE  MDHNKL  (Z ,H1 , H2 , H1PRME , H2PRME) 

G  COMPUTE  MODIFIED  HANKEL  FUNCTIONS  OF  ORDER  ONE  THIRD 

IMPLICIT  REAL  *8  (A-H.O-Z)  2 

COMPLEX*16  zsqrt.zexp 
REAL* 8  zabs 

COMPLEX* 16  Z , I , HI , H2 , H1PRME , H2PRME , ZPOWER , TERM1 , TERM2 ,  3 

$  TERM3 , ZTERM , TERM , SUM1 , SUM2 , SUM3 , SUM4 , SQRTZB ,  4 

$  EXP1 , EXP2 ,  EXP3 , EXP4 , EXP5 , GM2F , GPMFP , MPOWER , BETA , RTZ ,  5 

$  CONST1 , CONST2 , CONST3 , CONST4  6 

DIMENSION  A(23) ,  B(23) ,  C(23) ,  D(23),  CAP (14)  7 

DATA  A/  8 

$  9.30436716930000D-01.3.10145572309700D  01 , 2 . 06763714873160D  02,  9 

$  5 . 74343652425450D  02 , 8 . 70217655190080D  02 , 8 . 28778719228640D  02,  10 

$  5 .4168543 7404 340D  02,2.579454463830200  02,9.345849506631000  01,  11 

$  2 . 66263518707400D  01,6.121000430056000  00 , 1 . 15928038448000D  00,  12 

$  1 . 84012759441000D-01 , 2 . 48330309640000D-02 , 2 . 88420801000000D-03 ,  13 

$  2 . 91334142000000D-04, 2 . 582749 50000000D- 05 , 2 . 02568600000000D-06  ,  14 

$  1 . 41557000000000D-07 , 8 . 87000000000000D-09 , 5 . 01000000000000D-10 ,  15 

$  2 . 60000000000000D- 11 , 1 . 00000000000000D- 12/  16 

DATA  B/  17 

$  6 . 78298725140C00D-01 , 1 . 13049787524000D  01 , 5 . 38332321543100D  01,  21 

$  1 . 1962940478 73 50D  02 , 1 . 53371031778650D  02 , 1 . 27809193148880D  02,  19 

$  7 . 47422 182 15 720QD  01 , 3 . 23559386215200D  01 , 1 . 07853128738400D  01,  18 

$  2 . 85325737403000D  00 , 6 . 13603736351000D-01 , 1 . 09376780098000D-01 ,  23 

$  1 . 64229399550000D-02 , 2 . 10550512200000D-03 , 2 . 33167788000000D-04 ,  22 

$  2 . 25282890000000D-05 , 1 . 91567100000000D-06 , 1 . 44470000000000D-07 ,  20 

$  9 . 72900000000000D-09 , 5 . 89000000000000D- 10 , 3 . 20000000000000D- 11 ,  24 

$  2 . OOOOOOOOOOOUoOD-12 , 0 . OOOOOOOOOOOOOOD  00/  25 

DATA  C/  26 

$  4 . 65218358460000D-01 , 6 . 202911446190G0D  00 , 2 . 58454643591500D  01,  27 

$  5 . 22130593114000D  01 , 6 . 21584039421500D  01,4.875168936639000  01,  28 

$  2 . 70842 7 18 702 7 00D  01 , 1 . 12150194079600D  01 , 3 . 59455750255000D  00,  29 

$  9 . 18 1500645 10000D- 01 . 1.9128126 34 39 000D-01, 3. 31222966990000D-02,  30 

$  4 . 8424410 3 800C00D- 03 , 6 . 05683682000000D-04 , 6 . 55501820000000D-05  ,  31 

$  6 . 19859900000000D-06 , 5 . 16550000000000D-07 , 3 . 82200000000000D-08 ,  32 

$  2 . 52800000000000D-09 , 1 . 50000000000000D- 10 , 8 . OOOOOOOOOOOOOOD- 12 ,  33 

$  0. OOOOOOOOOOOOOOD  00,0. OOOOOOOOOOOOOOD  00/  34 

DATA  D/  35 

$  6 . 7829872 5140000D- 01 , 4 . 52199150096200D  01 , 3 . 76832625080150D  02,  36 

y  1. 19629404787350D  03,1.993823413122500  03 , 2 . 04494709038206D  03,  37 

$  1.420 102 14609865D  03 , 7 . 11830649673510D  02 , 2 . 69632821846030D  02,  38 

$  7. 989120 64729000D  01 , 1 . 90217158268800D  01 , 3 . 71881052333900D  00,  39 

$  6 . 07648778323000D-01 , 8 .42202048960000D-02 , 1 . 00262148690000D-02 ,  40 

$  1 .03630127800000D-03 , 9 . 38678690000000D-05 , 7 . 51243500000000D-06 ,  41 

$  5 . 35074000000000D-07 , 3 .41350000000000D-08 , 1 . 96200000000000D-09 ,  42 

$  1.0200C000000000D- 10, 5. OOOOOOOOOOOOOOD- 12/  43 

DATA  CAP/  44 

$  1.04166666666667D-01,8.35503472222222D-02,1.28226574556327D-01,  45 

$  2. 9184902 6464 140D-01, 8. 8162726744 3758D-01, 3. 32140828186 277D  00,  46 

$  1.49957629868626D  01 , 7 . 89230130115870D  01 , 4 . 74451538868000D  02,  47 

$  3 . 20749009100000D  03 , 2 . 40865496000000D  04 , 1 . 98923120000000D  05,  48 

$  1 . 79190200000000D  06 , 1 . 74843770000000D  07/  49 

C  50 

DATA  1/(0. DO, 1. DO)/ 

DATA  ROOT3/1. 73205080756888D  00/  52 

DATA  ALPHA/8 . 53667218838951D-01/  53 

DATA  CONSTl/(  2 . 58819045102522D-01 , -9 . 65925826289067D-01)/  54 

DATA  CONST2/(  2 . 58819045102522D-01 ,  9 . 65925826289067D-01)/  55 
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DATA  C0NST3/( -9.65925826289067D-01,  2 . 58819045102522D-01)/  56 

DATA  C0NST4/( - 9 . 65925826289067D-01 , -2 . 58819045102 5 22D- 01) /  57 

C  58 

ZPOWER-l.O  59 

SUM3-0.0  60 

SUM4-0.0  61 

ZMAG-zabs (Z)  62 

I F ( ZMAG  .GT.  4.2)  GO  TO  70  63 

IF(ZMAG  .GE.  3.2)  GO  TO  10  64 

N-12  65 

GO  TO  30  66 

10  I F ( ZMAG  .GE.  4.1)  GO  TO  20  67 

N— 15  68 

GO  TO  30  69 

20  N-23  70 

30  SUM1-0.  71 

SUM2-0.  72 

ZTERM=-Z**3/200.0  73 

DO  50  M-1,N  74 

SUM1-SUM1+A(M)*ZP0WER  75 

SUM2-SUM2+B (M) *ZPOWER  76 

SUM3=-SUM3+C(M)*ZPOWER  77 

SUM4=*SUM4+D(M)  *ZPOWER  78 

ZPOWER-ZPOWER*ZTERM  79 

IF(zabs (ZPOWER)  .LE.  1.0D-30)  GO  TO  60  80 

50  CONTINUE  81 

60  GM2F~I*(Z*SUM2-2.*SUMl)/R00T3  82 

GPMFP“I*(SUM4+2 . *Z*Z*SUM3)/ROOT3  83 

H1-Z*SUM2+GM2F  84 

H2=-H1  - 2 . 0*GM2F  85 

H1PRME-SUM4+GPMFP  86 

H2PRME-H1PRME- 2 . 0*GPMFP  8/ 

RETURN  88 

C  89 

70  SUM1-1.0  90 

SUM2-1.0  91 

RTZ-zsqrt(Z)  92 

SQRTZB-RTZ*Z  9j 

ZTERM-I/SQRTZB  94 

MPOWER-l.O  95 

TERM--1.5/Z  96 

DO  80  M-1,14  Q 7 

ZPOWER-ZPOWER*ZTERM  98 

MPOWER-MPOWER*(-ZTERM)  99 

TERMl^CAP (M)*ZPOWER  100 

TERM2-CAP(M)*MP0WER  101 

SUM1=SUM1+TERM1  102 

SUM2-SUM2+TERM2  103 

SUM3-SUM3+M*TERM1  104 

SUM4-SUM4+M*TERM2  105 

80  CONTINUE  106 

SUM3-SUM3*TERM  107 

SUM4-SUM4*TERM  108 

EXPl”Zexp(2 . *I*SQRTZB/3  .  )  109 

EXP2=EXP1*C0NST1  110 

EXP3-CONST2/EXP1  111 

EXP4=C0NST3*EXP1  112 

EXP5=C0NST4/EXP1  113 
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BETA-ALPHA/zsqrt(RTZ)  114 

ZREAL-Z  115 

ZIMAG—  I*Z  116 

IF  (ZREAL.GE. 0.0. OR. ZIMAG. GE.O.O)GO  TO  90  117 

HI— BETA* (EXP2*SUM2+EXP5*SUM1 )  118 

H1PRME-BETA*(EXP2*(SUM2*( -0 . 25/Z+I*RTZ)+SUM4)+EXP5*(SUMl*( -0 . 25/Z  119 

$  - I*RTZ ) +SUM3 ) )  120 

GO  TO  110  121 

90  H1-BETA*EXP2*SUM2  122 

H1PRME-BETA*EXP2* ( SUM2* (-0.25 /Z+I*RTZ) +SUM4 )  123 

110  IF  (ZREAL.GE. 0.0. OR. ZIMAG. LT.O.O)GO  TO  120  124 

H2— BETA* ( EXP3*SUM1+EXP4*SUM2 )  125 

H2PRME-BETA* ( EXP 3* ( SUM1* ( - 0 . 2  5/Z - I*RTZ ) +SUM3 ) +EXP4* ( SUM2* ( - 0 . 2  5/Z  126 

$  +I*RTZ)+SUM4) )  127 

RETURN  128 

120  H2-BETA*EXP3*SUM1  129 

H  2  PRME-B  ETA*  EXP  3  * ( SUM1* ( - 0 . 2  5/Z - I*RTZ ) +SUM3 )  130 

RETURN  131 

END  132 


SUBROUTINE  MAGANG(ARG , MAG .ANGLE) 

IMPLICIT  REAL  *8(A-H,0-Z) 

REAL* 8  DSQRT.dacos 
COMPLEX*16  ARG.IM 
REAL*8  MAG 

DATA  RDTDEG/5 . 729577951D+01/ 
data  IM/(0 . DO , 1 . DO) / 
c 

ENTRY  MGNGLE(ARG, MAG .ANGLE) 

ARGRAL  -  ARG 
ARGMAG  -  - IM*ARG 

MAG  -  DSQRT (ARGRAL* ARGRAL  +  ARGMAG*ARGMAG ) 

IF(MAG  .EQ.  0.0)  GO  TO  10 
COSQ  -  ARGRAL/MAG 

IF  (COSQ  .LT.  -1.0. AND.  COSQ  .GT.  -1.01)  COSQ  -  -1.0 
IF(COSQ  .GT.  1.0. AND.  COSQ  .LT.  1.01)  COSQ  -  1.0 
5  ANGLE  -  dacos ( COSQ) *RDTDEG 

IF (ARGMAG  .LT.  0.0)  ANGLE  -  360.0  -  ANGLE 
RETURN 

10  COSQ  -  0.0 
GO  TO  5 
END 

SUBROUTINE  CLIN  EQ  (A,  B,  X,  N.  N  DIM,  I FLAG ,  ERR) 

C 

C  CLIN  EQ  USES  L-U  DECOMPOSITION  TO 
C  FIND  THE  TRIANGULAR  MATRICES  L,  U 
C  SUCH  THAT  L  *  U  -  A.  L  AND  U  ARE 

C  STORED  IN  A.  THIS  FORM  IS  USED  WITH 

C  BACK- SUBSTITUTION  TO  FIND  THE  SOLN 

C  X  OF  A*X— L*U*X-B. 

C  N  IS  THE  NUMBER  OF  EQUATIONS  AND 
C  N  DIM  IS  THE  DIMENSION  OF  ALL  ARRAYS 
C  IN  THE  PARAMETER  LIST. 

C 

C  IF  I FLAG  -  0,  L,  U,  AND  X  ARE 
C  COMPUTED. 

C  IF  I FLAG  IS  NON-ZERO,  IT  IS  ASSUMED 
C  THAT  L  AND  U  HAVE  BEEN  COMPUTED  IN 
C  A  PREVIOUS  CALL  AND  ARE  STILL  STORED 
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C  IN  A.  THUS  ONLY  X  IS  COMPUTED. 
C  ERR  IS  THE  ESTIMATED  RELATIVE 
C  ERROR  OF  THE  SOLUTION  VECTOR. 

C 


COMPLEX*16  A,  B,  X,  T 

00002300 

I NT EGER* 2  IROW 

DIMENSION  A(N  DIM,  N  DIM), 

00002500 

$  B(N  DIM) ,  X(N  DIM) 

00002600 

DIMENSION  IROW(50) ,  Q(50) 

00002700 

DATA  EPS  /I . 0E- 15/ 

00002800 

IF  (N.GT.50)  GO  TO  900 

00002900 

00003000 

00003100 

IF  (IFLAG.NE.O)  GO  TO  600 

00003200 

DO  050  I  -  I.N 

00003300 

Q(I)  -  0.0 

00003400 

DO  040  J  -  I,N 

00003500 

QQ  -  zabs  (A(I , J) ) 

040  IF  (Q( I) . LT.QQ)  Q(I)  -  QQ 

00003700 

IF  (Q(I) .EQ.0.0)  GO  TO  901 

00003800 

050  CONTINUE 

00003900 

ERR  -  EPS 

00004000 

PPIV  -  0.0 

00004100 

DO  100  I  -  1 ,N 

00004200 

100  IROW(I)  -  I 

00004300 

DO  500  L  -  1 , N 

00004400 

00004500 

PIVOT  -  0.0 

00004600 

K  -  L  -  1 

00004700 

DO  240  I  -  L,N 

00004800 

IF  (K.LT.l)  GO  TO  230 

00004900 

DO  220  J  -  1 , K 

00005000 

220  A(I , L)  -  A( I , L)  -  A(J , L)  *  A(I , J) 

00005100 

230  F  -  zabs  (A(I,L))  /  Q(I) 

IF  (PIVOT. GT.F)  GO  TO  240 

00005300 

PIVOT  -  F 

00005400 

NPIVOT  -  I 

00005500 

240  CONTINUE 

00005600 

IF  (PIVOT. EQ.0.0)  GO  TO  901 

00005700 

IF  (PPIV. LE. PIVOT)  GO  TO  250 

00005710 

ERR  -  ERR  *  PPIV  /  PIVOT 

00005720 

TF  (ERR. GE. 1.0)  GO  TO  901 

00005730 

250  PPIV  -  PIVOT 

00005740 

IF  (NPIVOT. EQ.L)  GO  TO  280 

00005750 

Q(NPIVOT)  -  Q(L) 

00005800 

J  -  IROW(L) 

00005900 

IROW(L)  -  IROW(NPIVOT) 

00006000 

IROW (NPIVOT)  -  J 

00006100 

DO  260  I  -  1,N 

00006700 

T  -  A(L, I) 

00006800 

A(L, I)  -  A(NPIVOT , I ) 

00006900 

A(NPIVOT,I)  =  T 

00007000 

260  CONTINUE 

00007100 

280  IF  (L.EQ.N)  GO  TO  500 

00007110 

T  -  (1.0D0.0.0D0)  /  A(L,L) 

00007200 

K  -  L  +  1 

00007300 

M  -  L  -  1 

00007400 

DO  450  I  -  K , N 

00007600 
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IF  (M.LT.l)  GO  TO  400 
DO  350  J  -  1,M 

350  A(L, I)  -  A(L, I)  -  A(L, J)  *  A(J,I) 
400  A(L, I )  -  T  *  A(L, I) 

450  CONTINUE 
500  CONTINUE 

IF  (ERR.GT.1.0E-5)  PRINT  998,  ERR 


600  DO  620  I  -  2 ,N 

620  X(I)  -  (0.0D0.0.0D0) 

J  -  IROW(l) 

X(l)  -  B(J)  /  A(1 , 1) 

DO  700  I  -  2,N 
J  -  IROW(I) 

K  -  I  -  1 
DO  650  L  -  l.K 

650  X(I)  -  X(I)  +  A(I,L)  *  X(L) 

X(I)  -  (B(J)  -  X(I))  /  A(I , I) 

700  CONTINUE 
K  -  N  -  1 
DO  800  I  -  l.K 
J  -  N  -  I 
M  -  J  +  1 
DO  800  L  -  M,N 
X(J)  -  X(J )  -  X(L)  *  A(J , L) 

800  CONTINUE 
RETURN 

900  PRINT  999 
ERR  -  1.0 
RETURN 

901  PRINT  997 
ERR  -  1.0 
RETURN 

997  FORMAT  ( ' 1 ERROR  IN  CLIN  EQ,  MATRIX  IS  SINGULAR') 

998  FORMAT  ('  CAUTION-', 

$  '  CLIN  EQ  HAS  DECOMPOSED  AN  ILL-CONDITIONED  MATRIX.',/, 

$  '  RESULTS  WILL  HAVE  RELATIVE  ERROR  -'.E11.2) 

999  FORMAT  ( ' 1 ERROR  IN  CLIN  EQ,  MATRIX  SIZE  GREATER  THAN  50') 
END 
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SUBROUTINE  MCFLD 
parameter  (maxinds-30) 

IMPLICIT  REAL  *8<A-H,0-Z) 

C  COMPUTE  FIELDS  FROM  XVAL  MIN  TO  XVAL  MAX  FOR  TWO  XMTR-RCVR  DISTANCES 
C  AT  DELTAX  INTERVALS 
C 

COMMON/TERM/NT , NTR 

COMMON/MCINPT/THETA(25 .maxmds) , FOFR ( 2  5 , maxmds ) , 

$  XTRA(3, 3,25, maxmds) ,TOPHT(25) , 

$  XVAL(25) , FREQ, RHOMAX , RHOMIN , DELRHO , DELTAX ,epsr (25) , 

$  SIGMA(25) ,NRSLAB,NRMODE(25) ,NTMAX 

COMMON/MCSTOR/A(25 , maxmds .maxmds) , S (25 .maxmds) , C(25 .maxmds) , 

$  NTHSQ (25), KVRAOT , KVRATT , 

$  AVRKOT , AVRKTT .CONST, OMEGA , WAVENO 

character*68  idplot 

COMMON/MCPLOT /R( 400 ) ,DBy(400) ,DBx(400) ,ANGchi(400) , 

$  IDPLOT, ISUB 

COMMON/HTGN/F (4 , 25 .maxmds ,2) 

COMMON/HGINPT/GAMMA , PHI , ZT , ZR , S INGAM , COSGAM , S INPHI , COSPHI 

common/p 1 tflg/iplflg 

REAL*4  R, DBy , dbx , angchi 

real*8  num 

REAL* 8  KVRAOT , KVRATT 

complex*16  tahy , tahx 

COMPLEX* 16  SOLNA (maxmds, 3) , THETA , A , S , C , XTRA , TB , FOFR , F, NTHSQ , IM 
DATA  ERAD/6 . 370D3/ 
data  IM/ (0 . 0D0 , 1 . 0D0) / 

C 

C 

MP— 10 

RHO  -  RHO  MIN 

DO  1  LL=-2  ,  NRSLAB 

IF(XVAL(LL) -RHO  .GE.  0.)  GO  TO  2 

1  CONTINUE 
M-nrslab 
GO  TO  3 

2  M  -  LL-1 

3  CONTINUE 

IF (M  .EQ.  MP)  GO  TO  720 

DO  710  L-1,3 

DO  710  J  -  l.NRMODE(m) 

SOLN  A(J,L)  -  (0.0, 0.0) 

DO  710  K  -  1 ,NRMODE(ntr> 

IF(M  .NE.  NTR)  GO  TO  35 

if(l  .eq.  1)  then 

SOLN  A(J,L)  -  SOLN  A(J,L) 

$  +A(M, J ,K)*(XTRA(1 ,L,NTR,K)*F(1 ,NTR,K, 1)* 

$C0SGAM  +XTRA ( 2 , L , NTR , K) *F( 2 , NTR , K , 1 ) *S INGAM* 

$COSPHI  +XTRA( 3 , L , NTR , K)*F( 3 , NTR , K , 1 )*SINGAM* 

$SINPHI)/( -s(ntr.k)) 
else 

SOLN  A( J , L)  =  SOLN  A(J,L) 

$  +A(M, J ,K)*(XTRA(1 ,L,NTR,K)*F(1 , NTR.K, 1)* 

$COSGAM  +XTRA ( 2 , L , NTR , K) *  F ( 2 , NTR , K , 1 ) *S I NGAM* 

$COS  PH I  +XTRA ( 3 , L , NTR , K ) *F ( 3 , NTR , K , 1 ) *S INGAM* 

$SINPHI ) 
endif 
GO  TO  710 
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35  continue 

if(l  .eq.  1)  then 

SOLN  A(J,L)  -  SOLN  A(J,L) 

$  +A(M, J ,K)*(XTRA(1 , L,NTR,K)*F(1 ,NTR, K , 1)* 

$COSGAM  +XTRA ( 2 , L , NTR , K) *F ( 2 , NTR , K , 1 ) *S INGAM* 

$COSPHI  +XTRA(3,L,NTR,K)*F(3,NTR!K,1)*SINGAM* 

$SINPHI) 

$  *zexp( - IM*WAVENO*S ( NTR , K)*XVAL(NTR+ l))/(-s(ntr,k)) 
else 

SOLN  A(J,L)  -  SOLN  A(J,L) 

$  +A(M,J,K)*(XTRA(1,L,NTR,K)*F(1.NTR,K,1)* 

$COSGAM  +XTRA ( 2 , L , NTR , K) *F< 2 , NTR , K , I ) *S INGAM* 

$COSPHI  +XTRA( 3 , L, NTR , K)*F( 3 , NTR , K , 1)*S INGAM* 

$SINPHI ) 

$  *zexp( - IM*WAVENO*S (NTR , K) *XVAL(NTR+1) ) 
endif 

710  CONTINUE 
C 

720  CONTINUE 

tahy  -  (0.0, 0.0) 
tahx  -  (0.0,0 .0) 

DO  J-l,NRMODE(m) 

IF(M  .eq.  NTR)  then 

TB  -zexp ( - TM*WAVENO*S (M  ,J)*RH0) 
else 

TB  -zexp ( IM*WAVENO*S (M  ,J)*(XVAL(M  )  -  RHO)) 
endif 

tahy=tahy+SOLNA(J , 1)*TB*F(1 ,M, J ,2) 
tahx=tahx+SOLNA(J , 3)*TB*F(4,M, J , 2) 
enddo 

tahy-tahy*CONST/DSQRT (DSIN(RHO/ERAD) ) 

tahx-tahx*CONST/DSQRT (DS IN (RHO/ERAD) ) 

tahyr-tahy 

tahyi--im*tahy 

tahxr-tahx 

tahxi-- ira*tahx 

phshy-datan2 ( tahyi , tahyr ) 

phshx=-datan2  ( tahxi ,  tahxr ) 

amphy=-dsq  r  t  ( tahy  r  *  tahy  r+  tahy  i  *  tahy  i ) 

amphx-ds  q  r  t ( tahxr *  tahxr+  tahx i*  tahx i ) 

beta*=datan2  (amphy ,  amphx) 

de 1 ta-phshy - phshx 

num=ds in( 2 .0*beta)*dcos(delta) 

den— dcos (2 .0*beta) 

psi-0 . 5*datan2 (num , den) *57 .295779 

if(psi  .It.  0.0)  psi=  :>si+180 . 0 

chi-90 . 0-psi 

amp y- 10 . 0*dlogl0(dmaxl(l . Od- 30 , amphy* amphy )*1 . 0dl2) 
ampx-10 . 0*dlogl0 (dmaxl ( 1 . Od- 30 , amphx* amphx) *1 . 0dl2 ) 

R(NT)  -  XVAL( 2 ) 
dby(nt)  -  ampy 
dbx(nt)  -  ampx 
angchi(nt)  =  chi 
MP  -  M 

IF (NT  .NE.  NTMAX)  RETURN 
PGAMMA  -  GAMMA/1 .7453290-2 
PPHI  -  PHI/1. 745329D-2 
PRINT  92 7, PGAMMA, PPHI, ZT.ZR 
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927  FORMAT (//,*  GAMMA(DEG)-' , F6 . 1 , '  PHI (DEG)-' , F6 . 1 , '  ZT(KM)- 
$'  ZR(KM)-' ,F7.2) 
print  929 

929  FORMAT(/,13X, 'RHO-KM' ,3X, 'MAG(HY) -DB' ,3X, 'MAG(HX) -DB' ,3x, 

$  'DOA  ERROR -DEG' ) 

DO  JJ-l.NTMAX 

PRINT  908 ,R(JJ) , dby(j j ) , dbx( j j ) , angchi( j j ) 

908  FORMAT (9X , F10 . 2 , 2x , F10 . 5 , 3x , F10 . 5 , 3x , F10 . 5 ) 
enddo 

write(*, '(/)') 

IF  (iplflg. EQ.O)  RETURN 

CALL  MCPLTS 

RETURN 

END 


,  F7 . 2 , 
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SUBROUTINE  MCFLD2 
parameter  (maxmds-30) 

IMPLICIT  REAL  *8(A-H,0-Z) 

C  COMPUTE  FIELDS  FROM  RHO  MIN  TO  RHO  MAX 
C  AT  DEL  RHO  INTERVALS. 

C 

COMMON/TERM/NT , NTR 

COMMON/MCINPT/THETA(25 .maxmds) , FOFR(25 .maxmds) , 

$  XTRA(3 , 3 , 25 .maxmds) ,TOPHT(25) , 

$  XVAL( 25 ) , FREQ , RHOMAX , RHOMIN , DELRHO , DELTAX , epsr ( 25 ) , 

$  SIGMA(25) ,NRSLAB,NRMODE(25) .NTMAX 

COMMON/MCSTOR/A( 25 .maxmds .maxmds) , S(25 .maxmds) , C (25 .maxmds) , 

$  NTHSQ(25) , KVRAOT , KVRATT , 

$  AVRKOT.AVRKTT, CONST, OMEGA. UAVENO 

character*68  idplot 

COMMON/MCPLOT/R(400) ,DBy(400) ,DBx(400) ,ANGchi(400) . 

$  IDPLOT, ISUB 

COMMON/HTGN/F ( 4 , 2 5 , maxmds , 2 ) 

COMMON /HGINPT/GAMMA , PHI , ZT , ZR , SINGAM , COSGAM , SINPHI , COSPHI 
common/pL tflg/iplf lg 
REAL*4  R , DBy , dbx , angchi 
COMPLEXES  SOLNA ( maxmds , 3) 

COMPLEX*lb  THETA , A, S , C , XTRA , TB , IM , FOFR , F , NTHSQ 

complex*16  tahy.tahx 

REAL*8  KVRAOT , KVRATT 

real*8  num 

DATA  ERAD/6 . 370D3/ 

data  IM/(0 . DO , 1 . DO)/ 

C 

C 

ISUB  -  1 
DBMAX  -  -1000.0 
RHO  -  RHO  MIN 
M  -  NTR 
X  -  RHO  -  1.0 
600  CONTINUE 

700  IF  (RHO.LE.X)  GO  TO  720 
DO  710  L-1,3 
DO  710  J  -  1 . NRMODE(m) 

SOLN  A(J,L)  -  (0.0, 0.0) 

DO  710  K  -  1 .NRMODE(ntr) 

IF(M  . NE .  NTR)  GO  TO  35 

if ( 1  .eq.  1)  then 

SOLN  A(J,L)  -  SOLN  A(J,L) 

$  +A(M , J , K)*(XTRA(1 , L, NTR ,K)*F(1 , NTR, K , 1)* 

§COSGAM  +XTRA( 2 , L , NTR , K) *F( 2 , NTR , K , 1)*S INGAM* 

$COSPHI  +XTRA(3,L,NTR,K)*F(3,NTR,K,1)*SINGAM* 

$ SINPHI)/ (-s(ntr.k)) 
else 

SOLN  A( J , L)  -  SOLN  A(J,L) 

$  +A(M,J,K)*(XTRA(1,L,NTR,K)*F(1,NTR,K,1)* 

$COSGAM  +XTRA ( 2 , L , NTR , K) *F( 2 , NTR , K , 1)*SINGAM* 

$COSPHI  +XTRA( 3 , L.NTR, K)*F(3 , NTR , K, 1)*S INGAM* 

$SINPHI ) 
endif 
GO  TO  710 
35  continue 

if(l  .eq.  1)  then 
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SOLN  A( J , L)  -  SOLN  A(J,L) 

$  +A(M, J ,K)*(XTRA(1 , L,NTR,K)*F(1 , NTR,K, 1)* 

$COSGAK  +XTRA(2,L,NTR,K)*F(2,NTR,K,1)*SINGAM* 

$COSPHI  +XTRA( 3 , L , NTR , K) *F ( 3 , NTR , K , 1 ) *SINGAM* 

$SINPHI) 

$  *zexp ( - IM*WAVENO*S (NTR , K)*XVAL(NTR+1) )/ ( - s (ntr , k) ) 
else 

SOLN  A(J.L)  -  SOLN  A(J,L) 

$  +A(M, J ,K)*(XTRA(1 ,L,NTR,K)*F(1 ,NTR,K, 1)* 

$COSGAM  +XTRA (2 , L, NTR ,K)*F(2 , NTR , K , 1 ) *S INGAM* 

$COSPHI  +XTRA( 3 , L, NTR, K) *F(  3 ,NTR, K , 1)*S INGAM* 

$S INPHI ) 

$  *zexp ( - IM*WAVENO*S (NTR , K)*XVAL(NTR+1) ) 
endif 

710  CONTINUE 
M  -  M  +  1 
X  -  XVAL(M) 

if(m  .gt.  nrslab)  x— 1.0e6 
GO  TO  700 
C 

720  CONTINUE 

tahy  -  (0 . 0 , 0 . 0) 
tahx  -  (0 . 0 , 0 . 0) 

DO  J-l.NRMODE(m-l) 

IF (M- 1  .eq.  NTR)  then 

TB  -zexp(-IM*WAVENO*S(M-l,J)*RHO) 
else 

TB  -zexp(IM*WAVENO*S(M-l,J)*(XVAL(M-l)  -  RHO)) 
endif 

tahy  -  tahy+SOLN  A(J , 1)*TB*F(1 ,M-1 , J , 2) 
tahx  -  tahx+SOLN  A(J , 3)*TB*F(4 ,M- 1 , J , 2) 
enddo 

tahy  -  tahy*CONST/DSQRT(DSIN(RHO/ERAD) ) 

tahx  -  tahx*CONST/DSQRT(DSIN(RHO/ERAD) ) 

tahyr-tahy 

tahyi-- im*tahy 

tahxr=tahx 

tahxi=- im*tahx 

phshy“datan2 ( tahyi , tahyr ) 

phshx=datan2 ( tahxi , tahxr ) 

amphy=dsqrt ( tahyr*tahyr+tahyi*tahyi) 

amphx=dsqrt ( tahxr*tahxr+tahxi*tahxi) 

be  ta=datan2 ( amphy , amphx ) 

de 1 ta=phshy - phshx 

num=dsin(2 . 0*beta)*dcos (delta) 

den-dcos(2 .0*beta) 

psi=*0 . 5*datan2 (num,  den)*57. 295779 

if(psi  .It.  0.0)  psi«psi+180 . 0 

chi-90 . 0-psi 

ampy-10 . 0*dlogl0(dmaxl(l . Od- 30 , amphy*amphy)*l . 0dl2) 
ampx=10 . 0*dlogl0(dmaxl(l .  0d-30  ,  amphx* amphx )*1 . 0dl2) 

R(ISUB)  -  RHO 
DBy(ISUB)  -  ampy 
DBx(ISUB)  =  ampx 
angchi ( isub)-chi 
RHO  -  RHO  +  DEL  RHO 
ISUB  -  ISUB+1 

IF  (RHO  .LE.  RHOMAX)  GO  TO  600 
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I SUB  -  ISUB-1 
PGAMMA  -  GAMMA/1. 745329D- 2 
PPHI  -  PHI/1. 745329D-2 
PRINT  92 7, PGAMMA, PPHI ,ZT,ZR 
PRINT  929 

927  FORMAT (// , 11X, 'GAMMA (DEG)-' , F6 . 1 , '  PHI (DEG)-' , F6 . 1 , '  ZT(KM)- 
$  F7 . 2 , '  ZR(KM)-' ,F7.2) 

929  FORMAT (/ , 13X , 'RHO-KM' ,3X, 'MAG(HY)-DB' ,3X, 'MAG(HX) -DB' ,3x, 

$  'DOA  ERROR -DEG') 

DO  J-l.ISUB 

PRINT  908, R(J) , DBy(J) ,DBx(J) ,angchi(J) 

908  FORMAT (9X ,F10.2,2x,F10.5,3x,Fl0.5,3x, F10 . 5 ) 
enddo 

write (* , '(/)') 

IF  (iplflg.EQ.O)  RETURN 

CALL  MCPLT2 

RETURN 

END 
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SUBROUTINE  MCPLTS 
parameter  (maxmds-30) 

C  MCPLTS  GENERATES  TWO  PLOTS  (FIELD  AMPLITUDE  IN  DB  ABOVE  A 
C  MICRO  VOLT  PER  METER  FOR  1  KW  RADIATED  POWER  VERSUS  TRANSMITTER - 
C  TERMINATOR  DISTANCE  FOR  TWO  RECEIVER  POSITIONS). 

COMMON/MCPLOT/R(400) ,DBy(400) ,DBx(400) ,ANGchi(400) , 

$  IDPLOT(17) , ISUB 

COMMON/XPLOT /XMIN , xmax , Xt ic , YMIN ,ymax , Ytic , SIZEX , SIZEY 
COMMON /HG INPT /GAMMA , PHI , ZT , ZR , SINGAM , COSGAM , SINPHI , COSPHI 
COMMON/MCINPT/THETA(25 .maxmds) , FOFR( 25 .maxmds) , 

$  XTRA( 3 ,3,25, maxmds ) , TOPHT (25 ) , 

$  XVAL( 2  5 ) , FREQ.RHOMAX, RHOMIN , DELRHO , DELTAX , epsr ( 25) . 

$  S I GMA (25) , NRSLAB , NRMODE (25) , NTMAX 

dimension  xl(2),yl(2) 

COMPLEX* 16  THETA, FOFR.XTRA 

REAL* 8  XVAL , FREQ , RHOMAX , RHOMIN , DELRHO , DELTAX , eps r , SIGMA , TOPHT 
REAL* 8  SINGAM, COSGAM, SINPHI, COSPHI .GAMMA, PHI , ZT , ZR 
LOGICAL  UP (400 ) , upl ( 2 ) 
c 

xinc-(xmax-xmin) /sizex 
y inc- (ymax- ymin) /sizey 
DO  J-1,2 

upl(J)-. FALSE, 
enddo 

DO  J-1,400 

UP(J)-. FALSE, 
enddo 

CALL  pltbgn 

CALL  plot(l. 0,1.0, -3) 

call  symbol (2 . 0 , -0 . 5 , .10, 3 5HTransmit ter -Terminator  Distance(Mm) 
$  0.0,35) 

call  symbol ( -0 . 9 , 0 . 9 , . 10 , 26H  Bearing  Error  (Deg)  ,90.0,26) 
call  symbol ( -0 . 7 , 0 . 9 , . 10 , 35HAmplitude  (dB  above  luv/m  for  1  kW) 
$  90.0,35) 

CALL  borar2(sizex,XM IN/1000, xmax/1000 , Xt i c/ 1000, 2. *Xtic/ 1000,1, 
$  s izey , YMIN ,ymax , Ytic , 2 . *Ytic , - 1) 

GAMMAD  -  GAMMA/1 . 745329D-2 
PHID  -  PHI/1 . 745329D-2 
xpos-0 . 0 
xstar t-xmin 
c 

c  draw  legend 

xl  ( 1 ) =xs  tar t 
xl(2)=xstart+xinc*0. 7 
y 1 ( 1 ) -yuax+y inc*0 . 3 
yl(2)=ymax+yinc*0. 3 

call  curve (xl ,yl , upl , 2 , xmin , ymin, xinc ,y inc , 2 ) 
call  symbol (xpos+0 . 7 , s izey+0 . 25 , . 1 , 6H  |HY|,0.,6) 
c 

xl(l)-xl(l)+l. 5*xinc 
xl(2)~xl(2)+l. 3*xinc 

call  curve (xl , y 1 , upl , 2 , xmin , ymin , xinc , y inc , 4 ) 
xpos-xpos+1 . 5 

call  symbo 1 (xpos+0 . 5 , s izey+0 . 25 ,. 1 , 6H  |HX|,0.,6) 
c 

xl(l)-xl(l)+l. 5*xinc 
xl(2)-xl(2)+l. 5*xinc 

call  curve (xl ,y 1 , up 1,2, xmin, ymin, xinc ,yinc , 5) 


xpos-A^os+1  5 

call  symbol (xpos+O . 5 , sizey+0 . 25 , . 1 , 14H  Bearing  Error, 0., 14) 
c 

c  draw  solid  line  at  0.0 

xl(l)-xmin 
xl(2)-xmax 
yKD-c  C 
yl(2) -u . 0 

call  curve  (xl  ,yl,upl,2,xmin,ymin,  (xmax- xmin) /sizex , 

S  (ymax-ymin)/sizey , 1) 

c 

CALL  CURVE (R , dby , up , NTMAX , xmin , ymin , xinc , yinc , 2 ) 

CALL  CURVE (R , dbx , up , NTMAX , xmin, ymin, xinc, yinc, 4) 

CALL  CURVE (R, angchi , up .NTMAX , xmin , ymin, xinc ,yinc , 5) 

CALL  symbol (0 . 0 , sizey+1 . 3 , .10, IDPLOT ,0.0,68) 

CALL  symbol (0 ., sizey+0 . 9 ,. 10 , 17HFreq  -  kHz, 0.0, 17) 

CALL  number (0 . 7 , sizey+0 . 9 , . 10 , sngl ( FREQ) ,0.0,3) 

CALL  symbol( 3 . 5 , s izey+0 . 9 , . 10 , 29HReceiver  Distance  -  km, 

$  0.0,29) 

CALL  number ( 5 . 5 , sizey+0 . 9 , . 10 , sngl (rhomin) ,0.0,1) 

CALL  symbol(C ., sizey+0 . 7 ,. 10 , 29HZt  -  km  Zr  -  km, 

$  0.0,29)' 

CALL  number (0 . 50 , sizey+0 . 7 , . 10 , sngl(zt) ,0.0,2) 

CALL  number (2 . 05 , sizey+0 .7, .10,sngl(zr) ,0.0,2) 

CALL  symbol ( 3 . 5 , s izey+. 7 , . 10 , 

§  35HGamma  -  deg  Phi  -  deg, 0.0, 35) 

CALL  number (4 . 30 , sizey+0 . 7 , . 10 , ganunad ,0.0,1) 

CALL  number (6 . 05 , sizey+0 . 7 , . 10 , phid, 0 .0,1) 

CALL  pi tend 

RETURN 

END 

SUBROUTINE  MCPLT2 
parameter  (maxmds-30) 

C  MCPLT2  GENERATES  ONE  PLOT  (FIELD  AMPLITUDE  IN  DB  ABOVE  A  MICRO  VOLT 
C  PER  METER  FOR  1  KW  RADIATED  POWER  VERSUS  DISTANCE  FROM  TRANSMITTER) . 
COMMON/MCPLOT/R(400) ,DBy(400) ,DBx(400) ,ANGchi(400) , 

$  IDPLOT ( 17 ) , I SUB 

COMMON/XPLOT/XMIN,xmax,Xtic , YMIN , ymax , Ytic , SIZEX , SIZEY 
COMMON/HGINPT/GAMMA , PHI , ZT , ZR , SINGAM , COSGAM , SINPHI , COSPHI 
COMMON/MCINPT/THETA(25 .maxmds) , FOFR(25 .maxmds) , 

$  XTRA (3,3,25, maxmds ) , TOPHT (25), 

$  XVAL(25) , FREQ, RHOMAX, RHOMIN, DELRHO , DELTAX , epsr ( 25 ) , 

$  SIGMA(25) , NRSLAB ,NRMODE(25) , NTMAX 

dimension  xl(2),yl(2) 

COMPLEX* 16  THETA, FOFR, XTRA 

REAL+8  XVAL , FREQ , RHOMAX , RHOMIN , DELRHO , DELTAX , epsr , SIGMA , TOPHT 
REAL+8  S INGAM, COSGAM, SINPHI .COSPHI, GAMMA, PHI, ZT.ZR 
LOGICAL  UP (400) , upl ( 2 ) 
c 

xinc- (xmax- xmin) /sizex 
yinc-(ymax-ymin)/sizey 
DO  J-1,2 

upl(J)-. FALSE, 
enddo 

DO  J-1,400 

UP(J)-. FALSE, 
enddo 

CALL  pltbgn 
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c 

c 


c 


c 


c 

c 


CALL  plot (1.0, 1.0, -3) 

call  symbol ( 1 . 5 , -0 . 5 , . 10 , 21HDis tance  (Megame ters )  ,  0 . 0 , 21 ) 
call  symbol( -0 . 9 , 1 . 5 , . 10 , 26H  Bearing  Error  (Deg)  ,90.0,26) 
call  symbol ( -0 . 7 , 1 . 5 , . 10 , 35HAmplitude  (dB  above  luv/ra  for  1  kW) , 

$  90.0,35) 

CALL  bordr2 ( sizex , XMIN/1000 , xmax/1000 ,Xtic/1000 , 2 . *Xtic/1000 , - 1 , 

$  sizey , YMIN ,ymax,Ytic , 2 . *Ytic , -1) 

GAMMAD  -  GAMMA/1. 745329D-2 
PHID  -  PHI/1. 745329D-2 
xpos-0 . 0 
xstart=xmin 

draw  legend 
xl(l)-xstart 
xl ( 2 )-xs tart+xinc*0 . 7 
yl(l)-yma>'  -yinc*0.3 
yl(2)-ymax+yinc*0 . 3 

call  curve  (xl  ,yl,upl,2,xmin,ymin,xinc,yinc,2) 
call  symbol (xpos+0 . 7 , sizey+0 . 25 ,. 1 , 6H  |HY|,0.,6) 

xl(l)-xl(l)+l . 5*xinc 
xl(2)-xl(2)+l . 3*xinc 

call  curve (xl , yl , upl , 2 , xmin , ymin , xinc , yinc ,4) 
xpos— xpos+1 . 5 

call  symbol(xpos+0 . 5 , sizey+0 . 25 , . 1 . 6H  |HX|,0.,6) 

xl(l)-xl(l)+l . 5*xinc 
xl(2)-xl(2)+l . 5*xinc 

call  curve (xl ,yl , upl , 2 , xmin , ymin , xinc , yinc , 5) 
xpos-xpos+1 . 5 

call  symbol (xpos+0. 5 , sizey+0. 25 ,. 1 , 14H  Bearing  Error, 0., 14) 

draw  solid  line  at  0.0 

xl(l)=xmin 

xl(2)=xmax 

yl(l)-0.0 

yl(2)-0.0 

call  curve (xl , yl ,upl , 2 , xmin, ymin, (xmax- xmin) /sizex, 

$  (ymax-ymin)/sizey , 1) 

CALL  CURVE (R , dby , up , I SUB , xmin , ymin , xinc , y inc , 2 ) 

CALL  CURVE (R , dbx , up , I  SUB , xmin , ymin , xinc , y inc , 4 ) 

CALL  CURVE (R , angchi , up , I SUB , xmin , ymin , xinc , yinc , 5 ) 

CALL  symbol (0 . 0 , s izey+1 . 5 , . 10 , IDPLOT ,0.0,68) 

CALL  symbol ,0 ., s izey+1 . 1 ,. 10 , 17HFreq  -  kHz, 0.0, 17) 

CALL  number (0 . 7 , s izey+1 . 1 , . 10 , sngl ( FREQ) ,0.0,3) 

CALL  symbol (0 ., sizey+0 . 9 ,. 10 , 29HZt  =  km  Zr  -  km, 

$  0.0,29) 

CALL  number (0 . 50 , sizey+0 . 9 , . 10 , sngl (zt) ,0.0,2) 

CALL  number (2 . 05 , sizey+0 . 9 , . 10 , sngl (zr ) ,0.0,2) 

CALL  symbol (0 ., sizey+0 . 7 ,. 10 , 33HGamma  -  deg  Phi  -  deg, 

$  0.0,33) 

CALL  number ( 0.80, sizey+0 . 7 , . 10 , gammad ,0.0,1) 

CALL  number ( 2 .45, sizey +0.7, . 10, ph id, 0.0,1) 

CALL  pltend 

RETURN 

END 
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subroutine  pltbgn 
character*l  answr 
logical  first, autopl 
dimension  ia(8) , ibuff (2) 
data  ia/82 ,79,57,48,73,87,73, 80/ 
c  ASCII  R090IWIP 

data  first/ . true ./ , autopl/ . false . / 
if(first)  then 

open(unit-52 , f ile-' /dev/plt7550a' ) 

print  *,'If  this  is  the  hp  7550  plotter  and  you  want  auto  feed,' 
print  *,'then  set  up  the  plotter,  load  a  sheet  and  answer  y : ' 
print  * , ' Do  you  want  auto  feed?' 
read  1 ,  answr 
1  format (al) 

if (answr  .eq.  'y'  .or.  answr  . eq .  'Y')  then 
autopl- . true . 
else 

autopl- . false . 
end  if 
end  if 

if ( .not. autopl  .or.  first)  then 

print  *,'Set  up  plotter,  enter  rotation  (y/n)  when  ready' 
read  1 ,  answr 
first-. false . 
end  if 

call  hpinit(2 , 0 , 0 , 0 , 52) 

if (answr  .eq.  'y'  .or.  answr  .eq.  'Y') 

$  call  buff (1 , ia.xbuff , 8) 
call  newpen(l) 
return 
end 

subroutine  plton 
dimension  ia(3) 
data  ia/27 , 46 , 89/ 
c  ASCII  esc  .  Y 

call  buff (1 , ia.xbuff , -3) 

call  newpen(l) 

return 

end 

subroutine  pltend 

call  newpen(O) 

entry  pltoff 

call  plot (0 . 0 , 0 . 0 , 999) 

return 

end 


i 
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c 


subroutine  bordr2(xlng,xmin,xmax,xticl ,xtic2 ,ndx, 
$  ylng,ymin,ymax,yticl ,ytic2 ,ndy) 


xscale-xlng/(xmax-xmin) 

yscale-ylng/(ytnax-ymin) 

if (xticl*xscale  .le.  0.  .or.  xtic2*xscale  . le .  0.)  go  to  999 
if (yticl*yscale  .le.  0.  .or.  ytic2*yscale  .le.  0.)  go  to  999 

if(iabs(ndx)  ,gt.  9)  then 
sx- . 15 

nx-ndx- (ndx/10)*10 
else 
sx- .  1 
nx=ndx 
end  if 
xo-. 5*sx 

if(iabs(ndy)  . gt .  9)  then 
sy=. 15 

ny-ndy- (ndy/10)*10 
else 
sy- .  1 
ny-ndy 
end  if 
yo-. 5*sy 
c 

xres-abs (xticl)/10 . 
yres-abs (yticl)/10 . 
c 

tl— . 05 
t2- . 10 
yval-ymin 
ytc2-yrain 
xp-0 . 
yp-0. 
go  to  115 

112  yval-yval+yticl 

if (abs(yval-ymax)  .le.  yres)  then 
call  plot(0. ,ylng, 2) 
if(abs(yval-ytc2)  . le .  yres)  then 
xln--sy*(3+ny) 
yln-ylng-yo 

if(abs(yval)  .gt.  yres  .or.  abs(yval)  . ge .  10.) 

$  xln-xln-sy*aint(alogl0(abs(yval) ) ) 

if(abs(yval)  .It.  yres)  yval-0 . 
if(yval  .It.  0.)  xln-xln-sy 
call  plot(xln,yln, 3) 
call  number (xln.yln, sy ,yval , 0 . ,ny) 
end  if 

call  plot(0 . ,ylng, 3) 
go  to  120 
end  if 

yp-(yval-ymin)*y scale 
call  plot(xp ,yp , 2) 

if (abs(yval-ytc2)  .gt.  yres)  go  to  118 
call  plot(t2 ,yp , 2) 

115  xln--sy*(3+ny) 
yln-yp -yo 

if(abs(yval)  .gt.  yres  .or.  abs(yval)  .ge.  10.) 
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$  xln-xln-sy*aint(aloglO (abs (yval) ) ) 

if (abs (yval)  .It.  yres)  yval-O. 
if(yval  .It.  0.)  xln-xln-sy 
call  plot(xln, yin, 3) 
call  number (xln.yln , sy ,yval ,0 . ,ny) 
ytc2-ytc2+ytic2 
go  to  119 

118  call  plot( tl ,yp , 2) 

119  call  plot(xp ,yp , 3) 

if (abs (yval -ymax)  ,gt.  yres)  go  to  112 
c 

120  yp-ylng 
tl-ylng- . 05 
t2-ylng- . 10 
xval-xmin 
xtc2~xmin+xtic2 

122  xval-xval+xticl 

if (abs(xval-xmax)  .gt.  xres)  go  to  123 
call  plot(xlng,ylng, 2) 

if (abs(xval-xtc2)  .le.  xres)  xtc2-xtc2+xtic2 
go  to  130 

123  xp-(xval-xmin)*xscale 
call  plot(xp ,yp , 2) 

if (abs(xval-xtc2)  .gt.  xres)  go  to  128 
call  plot(xp , t2 , 2) 
xtc2-xtc2+xtic2 
go  to  129 

128  call  plot (xp , tl , 2 ) 

129  call  plot(xp ,yp , 2) 

if (abs(xval-xmax)  .gt.  xres)  go  to  122 
c 

130  xp-xlng 
tl-xlng- . 05 
t2-xlng- . 10 
ytc2-ytc2 -ytic2 

if (abs (ytc2 -ymax)  . le .  yres)  go  to  135 

132  yval-yval-yticl 

if (abs (yval -ymin)  .gt.  yres)  go  to  133 
call  plot(xlng,0 . , 2) 
go  to  140 

133  yp-(yval-ymin)*yscale 
call  plot(xp ,yp , 2) 

if (abs(yval-ytc2)  .gt.  yres)  go  to  138 
call  plot(t2,yp,2) 

135  ytc2-ytc2-ytic2 
go  to  139 

138  call  plot(tl ,yp , 2) 

139  call  plot(xp , yp , 2) 

if (abs (yval -ymin)  .gt.  yres)  go  to  132 
c 

140  yp-0 . 
tl-.05 
t2- . 10 
yin-- 2 ,*sx 
xtc2-xtc2 -xtic2 

if (abs(xtc2-xmax)  . le .  xres)  go  to  145 
142  xval-xval-xticl 

if (abs (xval-xmin)  . le .  xres)  then 
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call  plot(0 . , 0 . ,2) 
if (abs (xval-xtc2)  . le .  xres)  then 
xln--xo*(2+nx) 

if(abs(xval)  .gt.  xres  .or.  abs(xval)  . ge .  10.) 

xln-xln-xo*aint(aloglO(abs (xval) ) ) 
if(abs(xval)  .It.  xres)  xval-0 . 
if (xval  .It.  0.)  xln-xln-xo 
call  plot(xln,yln, 3) 
call  number (xln.yln, sx , xval ,0 . ,nx) 
end  if 

call  plot(0 . , 0 . ,3) 
return 
end  if 

xp- (xval - xmin) *xscale 
call  plot(xp ,yp , 2) 

if (abs (xval -xtc2)  . gt .  xres)  go  to  148 
call  plot(xp , t2 , 2) 
xln-xp - xo* ( 2+nx ) 

if (abs (xval)  .gt.  xres  .or.  abs (xval)  . ge .  10.) 

i  xln-xln-xo*aint(aloglO(abs(xval) ) ) 
if(abs(xval)  .It.  xres)  xval-0. 
if (xval  .It.  0.)  xln-xln-xo 
call  plot(xln,yln, 3) 
call  number (xln, yin ,sx, xval, 0. ,nx) 
xtc2-xtc2-xtic2 
go  to  149 
call  plot(xp , tl , 2) 
call  plot(xp ,yp , 3) 

if (abs (xval -xmin)  . gt .  xres)  go  to  142 

print  100 ,xlng, xmin, xraax, x t icl ,xtic2,ndx, 

J  ylng,ymintyinax,yticl,ytic2,ndy 

format ('OError  in  B0RDR2 : ' / 

?  ' Oxlng ,  xmin,  xmax,  xticl,  xtic2,  ndx  -  ' , IpSell . 3 , i5/ 

}  ' Oylng ,  ymin,  ymax,  yticl,  ytic2,  ndy  -  ’ , lp5ell . 3 , i5) 

call  pltend 


subroutine  curve (x,y ,up , nrpts , xmin.ymin.xinc ,yinc , line) 
c 

c  x,y,up  must  be  dimensioned  at  least  nrpts 
c  xmin.ymin  are  x,y  origin  in  user  units 
c  xinc.yinc  are  x,y  scales  in  user  units  per  inch 
c 


c  line-1 

solid 

c  2 

long  dash 

c  3 

medium  dash 

c  4 

short  dash 

c  5 

dotted 

c  6 

short  4-  long  dash 

c  7 

short  4-  short  +  long  dash 

c 

logical  up,upl,up2 

dimension  ipen(8) , joc(7) ,x(nrpts) ,y(nrpts) ,up(nrpts) 
data  ipen/2 , 2 , 2 , 3 , 2 , 3 , 2 , 3/, joc/18 ,  11,  14,  23,  32,  41,  16/ 
data  delr/.l/ 
c 

if (nrpts  .le.  1)  go  to  99 
c 

if (line)  1,2,3 

1  kk-mod (line , 7)+7 
go  to  4 

2  kk-0 

go  to  4 

3  kk-mod (line , 7) 

4  kk-kk+1 
jo-joc(kk)/10 

j  c-j  oc (kk) - 10*j  o 
ip-ipen( jo) 
c 

j-0 
dr-O. 
rhol-0 . 
rho2-delr 

pxl-(x(l) -xmin) /xinc 
pyl-(y(l) -ymin)/yinc 
upl-up(l) 

if (.not.  upl)  then 
c 

c  go  to  first  position  with  pen  up 
call  plot(pxl , pyl , 3) 
if(kk  .eq.  6)  then 
px2-(x(2) -xmin)/xinc 
py2-(y(2) -ymin)/yinc 
delx-px2-pxl 
dely-py2 -pyl 

rho-sqrt(delx**2+dely**2) 

if(rho  .eq.  0.)  then 
dx  6-delx*.l 
dy  6-dely*.l 
else 

dx  6-delx*delr/rho* . 1 
dy  6-dely*delr/rho* . 1 
end  if 

call  plot(pxl4-dx6  ,pyl4dy6 , 2) 
end  if 
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end  if 


c 

do  40  i-2,nrpts 
px2-(x(i) -xmin) /xinc 
py2-(y(i) 'ymin)/yinc 
up2-up(i) 
if(up2)  then 
dr-0. 
rhol-0 . 
rho2-delr 
go  to  39 
end  if 
if(upl)  then 

c  pen  has  been  up,  prepare  to  lower  pen 
call  plot(px2 ,py2 , 3) 
go  to  39 
end  if 

if(kk  .eq.  2)  go  to  38 
delx-px2 -pxl 
dely~py2 -pyl 

rho-s  q  r  t ( de lx**2+de ly**  2 ) 
rhol-rhol+rho 

if(rho2  .gt.  rhol)  go  to  38 

de lx-de lx*de lr/rho 

de ly-de ly*de 1 r /rho 

dx  6-delx*.l 

dy  6-dely*.l 

if (dr  .eq.  0.)  go  to  20 

dx-delx*dr/delr 

dy— dely*dr/delr 

pxl-pxl+dx 

pyl-pyl+dy 

go  to  21 

20  if(rho2  .gt.  rhol)  go  to  38 

pxl-pxl+delx 
pyl-pyl+dely 

21  call  plot(pxl ,pyl, ip) 

if(kk  .eq.  6)  call  plot(pxl+dx6 ,pyl+dy6 , 2) 

j-j+1 

ip-ipen(jo+mod(j , jc)) 
rho2-rho2+delr 
go  to  20 

38  call  plot(px2 ,py2 , ip) 
dr-rho2-rhol 

39  pxl-px2 
pyl-py2 
up 1 -up 2 

40  continue 

99  return 

end 
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